AD-A218  460 


OTIC  FILE  COPY 


EPORT  DOCUMENTATION  PAGE 


.  CLASSIFICATION  AUTHORITY 

NOT  APPLICABLE 


2b.  DECLASSIFICATION /DOWNGRADING  SCHEDULE 

NOT  APPLICABLE 


A.  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 


6a.  NAME  OF  PERFORMING  ORGANIZATION 

Defense  Mapping  Agency 


■SiiTsh  vESti  TSfi 


6c.  ADDRESS  (O'ty,  State,  and  ZIP  Code) 

12100  Sunset  Hills  Road  Suite  200 
Reston,  Virginia  22090-3207 


I  lb.  RESTRICTIVE  MARKINGS 

I  NONE _ 

3.  DISTRIBUTION;  availability  OF  REPORT 

Approved  for  public  release, 
distribution  unlimited 


5  MONITORING  ORGANIZATION  REPORT  NUMB£R(S) 

SAME 


SOL  I  7a  NAME  OF  MONITORING  ORGANIZATION 


7b  ADDRESS  (Ofy,  State,  and  ZIP  C 


8a.  NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 


8c.  ADDRESS  (City,  State,  and  ZIP  Code) 


8b.  OFFICE  SYMBOL  9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUUR 
(If  applicable) 

NOT  APPLICABLE 


10.  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 

PROJECT 

TASK 

NOT  APPLICABLE 

ELEMENT  NO 

NO. 

NO. 

N/A 

N/A 

N/A 

WORK  UNIT 
ACCESSION  NO. 


11.  TITLE  (Include  Secunfy  Classification) 

Guidelines,  Operational  Procedures  and  Quality  Control  for  the  Estimation  of  Geodetic  Point 
Positions  from  GPS  Data  Collected  with  the  TI4100  Receiver 


12,  PERSONAL  AUTHOR(S) 

Malys,  Stephen  and  Jensen,  Peter  A. 


13a.  TYPE  OF  REPORT 

Technical 


13b.  TIME  COVERED 
FROM  NM  TO  N/A 


14.  DATE  OF  REPORT  {Year,  Month,  Day)  15.  PAGE  COUNT 

1989,  November  13 


16.  SUPPLEMENTARY  NOTATION 

Prepared  to  accorrpany  EMA-developed  software  known  as  STARPREP  and  GASP 


17,  COSATI  CODES 


FIcLD  GROUP  SUB-GROUP 


18  SUBJECT  TERMS  (Continue  on  reverse  if  necessary  and  identify  by  block  number)  , 

Geodesy  *  Satellite  Geodesy  .  f  >  i"'  ! 

GPS  y  WGS  B4  J 


19  ABSTRACT  {Continue  on  reverse  if  necessary  and  identify  by.fa/ock  number) 

"~j  Hjo  j 


The  Texas  Instrumeni.';  4100  (TI4100)  receiver  collects  data  from  the  US  Defense  DepartmenTs^AVSTAK  Global  Posinoning  System  of 
Satellites  (GPS).  The  Defense  Mapping  Agency  (DMA)  and  a  number  of  other  organizatrons  have  utilized  these  receivers  to  exploit  GPS  for 
a  variety  of  geodetic  applications.  This  report  provides  a  comprehensive  set  of  operabona]  procedures  and  guidelines  for  the  csbrrration  of 
geodetic  point  positions  (absolute  positions)  from  data  collected  with  the  TI4 100  receiver.  The  specific  algorithms  which  are  addressed  here 
have  been  developed  into  efficient,  portable  FORTRAN  77  programs  which  are  known  as  STARPREP  and  GASP.  The  name 
STARPREP  is  a  contr..vUon  for  GEOSTaR  PREProcessor  while  GASP  is  an  acronym  tor  Geodetic  Absolute  Sequential  Positioning 
program.  These  two  programs  are  designed  to  work  in  concert  to  exploit  GPS  data  for  geodetic  point  positioning  applications.  A  few  other 
utility  programs  are  also  addressed  which  perform  functions  such  as  data  transfer,  data  screening  and  data  analysis.  A  set  of  minimum 
observation  requirements  are  specified  such  that  a  geodetic  point  positioning  survey  can  be  planned  and  conducted  efficiently. 

This  report  does  not  address  the  use  of  GPS  data  which  has  been  subject  to  Selective  Availability  (SA)  or  Anti-Spoofing  (AS).  This 
report  will  be  updated  when  guidelines  and  standard  operating  procedures  have  been  established  for  such  data. 


20.  DISTRIBUTION /AVAILABILITY  OF  ABSTRACT 

□  UNCLASSIFIED/UNLIMITEO  Cjt  SAME  AS  RPT  D  OTIC  USERS 

21.  ABSTRACT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 

22a.  NAME  OF  RESPONSIBLE  INDIVIDUAL 

Stephen  Malys 

22b  TELEPHONE  f/nclucfe  Area  Code) 
(202)  227-2154 

22c  OFFICE  SYMBOL 
DMASC/EG 

DD  FORM  1473, 84  mar 


83  APR  edition  may  be  used  until  exhausted 
All  otner  editions  are  obsolete 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


CLEARED 

FOR  OPEN  PUBLICATION 

DEC  2  8  881'  12 

DIRECTORiTt  FOR  FREEDOM  OF  |•'FORMAT^OM 
ANOSEr  --  “  ivitr  fOASD-PA) 

OEPA  a  LF  LLFENSE 


GUIDELINES, 

OPERATIONAL  PROCEDURES,  AND 
QUALITY  CONTROL 

FOR 

THE  ESTIMATION  OF  GEODETIC  POINT 
POSITIONS  FROM  GPS  DATA  COLLECTED 
WITH  THE  TI4100  RECEIVER 


90  02  21  077 


895  466 


Table  of  Contents 


I 

i 

I 


page 

1 .  Introduction  1 

1.1  GPS  Geodetic  Point  Positioning  Applications  1 

1.2  Processing  Flow  Overview  3 

2.  Mission  Planning  and  Observation  Requirements  5 

3.  Field  Validation  8 

4.  Translation  of  GESAR  or  BEPP/CORE  data  to  the  FICA  format  9 

4.1  Overview  9 

4.2  Data  Translation:  Operating  Procedures  9 

4.3  FICA  Files  from  9-Track  Tapes  13 

5.  Screening  and  Editing  FICA  files  14 

6.  Preprocessing  (STARPREP)  16 

6.1  Overview  16 

6.2  Sample  STARPREP  Runstreams  16 

6.3  Broadcast  or  Precise  Ephemerides  and  Satellite  Clock  States  18 

6.4  Point  File  19 

6.5  Quality  Control  -  Analysis  of  STARPREP  Output  21 

7.  Point  Position  Estimation  (GASP)  22 

7.1  Overview  and  Sample  Runstream  22 

7.2  Algorithm  Description  23 

7.3  Pseudorange  Minus  Carrier  Beat  Fliase  Biases  24 

7.4  Reference  Satellite  Sequencing  24 

7.5  Rejected  Data  24 

7.6  Batch  Least  Squares  Versus  Sequential  Estimation  25 

7.7  V aiiance-Covariance  and  Correlation  Matrices  25 

7.8  Quality  Control  -  Analysis  of  GASP  Output  25 

8.  Precision  and  Accuracy  28 

8.1  Definitions  28 

8.2  Reporting  positioning  results  28 

9.  Summary  30 

Acknowledgements  30 

List  of  References  /  Bibliography  3 1 

Appendix  A.  Field  Logs  and  File  Management  Sheet  32 

Appendix  B.  Sample  STARPREP  Output  42 


Appendix  C.  Sample  GASP  Output 
Appendix  D.  Troubleshooting 


61 

84 


List  of  Figures 


Figure  1 .  Relative  and  Absolute  Positioning  2 

Figure  2.  Generalized  Data  Flow  Diagram  3 

Figure  3.  Mission  Planning  Aids  6 

Figure  4.  Broadcast  and  Precise  STARPREP  Runstreams  17 

Figure  5.  STARPREP  Runstreams  for  Excluding  a  Satellite  18 

Figure  6.  Sample  of  a  Portion  of  a  Point  File  19 

Figure  7.  Point  File  Format  20 

Figure  8.  Sample  GASP  Runstreams  22 

Figure  9.  Formation  of  GASP  Observables  23 

Figure  10.  Graphic  Example  of  Accuracy  and  Precision  28 


List  of  Tables 

Table  1 .  Minimum  Observation  Requirements  for  Geodetic  Point  Positioning  with  the  TI4100  7 
Table  2.  Data  Corrections  and  Their  Functional  Dependencies  16 


1.  INTRODUCTION 


The  Texas  Instruments  4100  (TI4100)  receiver  eoUects  data  from  the  US  Defense  Department’s 
NAVSTAR  Global  Positioning  System  of  Satellites  (GPS).  The  Defense  Mapping  Agency  (DI^)  and  a 
number  of  other  organizations  have  utilized  these  receivers  to  exploit  GPS  for  a  variety  of  geodetic 
applications.  This  report  provides  a  comprehensive  set  of  operational  procedures  and  guidelines  for  the 
estimation  of  geodetic  point  positions  (absolute  positions)  from  data  collected  with  the  114 1(X)  receiver. 
The  specific  algorithms  which  are  addressed  here  have  been  developed  into  efficient,  potable  FORTRAN  77 
programs  which  are  known  as  STARPREP  and  GASP.  The  name  STARPREP  is  a  contraction  for 
GEQSTAR  PREProcessor  while  GASP  is  an  acronym  for  Qeodetic  Absolute  Sequential  Eositioning 
program.  These  two  programs  are  designed  to  work  in  concert  to  exploit  GPS  data  for  geodetic  point 
positioning  applications.  A  few  other  utility  programs  are  also  addressed  which  perform  functions  such  as 
data  transfer,  data  screening  and  data  analysis.  A  set  of  minimum  observation  requirements  are  specified 
such  that  a  geodetic  point  positioning  survey  can  be  planned  and  conducted  efficiently. 

This  report  does  not  address  the  use  of  GI^  data  which  has  been  subject  to  Selective  Availability 
(SA)  or  Anti-Spoofing  (AS).  This  report  will  be  tqxlated  when  guidelines  and  standard  operating  procedures 
have  been  established  for  such  data. 

1.1  GPS  Geodetic  Point  Positioning  Applications 

For  the  last  few  decades,  DMA  has  exploited  the  Navy  Navigation  Satellite  System  (NNSS) 
(TRANSIT,  Doppler)  for  geodetic  point  positioning.  Over  6000  geodetic  point  positions  have  been 
estimated  in  a  global  reference  frame  and  have  been  effectively  utilized  for  activities  such  as  mapping 
control,  estimation  of  transformation  parameters  (regional  datums  to  World  Geodetic  Systems),  and  other 
MC&G  activities.  Since  the  NNSS  is  scheduled  to  terminate  operations  when  GPS  becomes  fully 
operational,  and  since  there  is  a  continuing  requirement  to  generate  geodetic  point  positions,  DMA  has 
developed  an  ability  to  utilize  GPS  data  for  the  estimation  of  individual  geodetic  point  positions. 

While  many  organizations  have  developed  algorithms  and  programs  for  relative  positioning,  only 
DMA  and  the  US  Naval  Surface  Warfare  Center  (NSWC)  have  developed  proven  methods  to  exploit  GPS 
for  point  positioning.  The  relative  positioning  (differential  positioning)  and  point  positioning  (absolute 
positioning)  concepts  are  shown  graphically  in  Figure  1.  Relative  positioning  can  be  described  as  the 
estimation  of  the  vector  (baseline)  connecting  a  known  station  to  an  unknown  station  while  point 
positioning  can  be  described  as  the  estimation  of  the  vector  connecting  the  origin  of  a  global  reference  frame 
with  an  unknown  station  anywhere  on  the  surface  of  the  Earth.  The  NSWC  GPS  point  positioning 
algorithm  has  been  described  by  Hermann  [1987]  while  the  DMA  algorithm  has  been  presented  in  Malys 
and  Ortiz  [1989]  and  Malys  and  Jensen  [1989]. 

Before  planning  a  tracking  campaign,  geodetic  users  of  GPS  data  must  decide  on  the  accuracy  and 
precision  goals  of  the  survey  as  well  as  the  logistics  requirements  of  the  personnel  and  eo’  ipmenL  In 
general,  a  relative  positioning  survey  requites  a  large  degree  of  coordination  among  field  party  members 
since  simultaneous  tracking  is  essential  for  the  success  of  all  baseline  estimation  algorithms.  In  contrast, 
point  positmning  surveys  can  be  conducted  with  little  or  no  coordination  among  field  parties.  Any  relative 
positioning  survey  requires  a  minimum  of  two  receivers  while  point  positioning  surveys  can  be  conducted 
with  a  single  receiver. 

While  most  relative  positioning  software  packages  offer  the  potential  a  centimeter-level  baselines, 
the  coordinates  assigned  to  the  unknown  end  of  a  baseline  are  limited  in  accuracy  by  the  accuracy  of  the 
’known”  (fixed)  end  of  the  baseline.  When  the  fixed  end  of  a  baseline  is  a  very  well-known  station  (a 
fiducial  site)  the  baseline  to  the  unknown  station  can  be  estimated  with  sub-centimeter  accuracy  if  special 
care  is  given  to  cycle-slip  repair  and  atmospheric  modeling.  If  on  the  other  hand,  the  station  at  the  end 
of  a  baseline  is  not  very  well-known,  the  baseline  will  suffer  scale  distortion  at  a  level  proportional  to  the 
error  in  the  assumed  fixed-site  coordinates  [Hilla,  personal  communication,  1988].  In  many  cases,  reliable 
error  estimates  for  the  coordinates  of  the  fixed  end  of  a  baseline  will  not  be  available. 

Since  GPS  point  positioning  is  generally  simpler  and  less  costly  in  terms  of  logistics,  many 
geodetic  GPS  users  will  find  it  desirable  to  design  and  conduct  point  positioning  surveys  instead  of  the 
more  complicated  relative  positioning  surveys.  This  is  particularly  true  when  the  positioning  results  will 
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be  used  as  as  basis  for  mapping  controL  Moreover,  when  GPS  positions  are  to  be  used  in  the  estimation  of 
transformation  parameters  between  the  World  Geodetic  System  1984  (WGS  84)  and  a  local  datum,  these 
GPS  positions  must  be  expressed  in  the  WGS  84  reference  frame.  Once  the  NNSS  is  removed  from 
service,  the  only  direct  method  of  obtaining  WGS  84  station  coordinates  will  be  by  estimating  point 
positions  with  GPS  data.  Any  other  methods  will  be  subject  to  a  complicated  sequence  of  error 
propagations  which  will  undoubtedly  be  difficult  to  perfonn  with  any  degree  of  confidence. 


ABSOLUTE  POSITIONING  MLAHWE  IPOSITinKDMN© 
(POINT  POSITIONING)  (MTOIRISNTniAIL  ]P(Q)Sinni©MNG) 


Figure  1.  Relative  and  Absolute  Positioning 
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1.2  Processing  Flow  Overview 

A  streamlined  process  has  been  developed  for  the  flow  of  tracking  data  from  the  field  to  the  DMA 
point  positioning  program  (GASP).  This  process  is  summarized  below  in  Figure  2. 


Figure  2.  Generalized  Data  Flow  Diagram 

The  flow  diagram  depicted  in  Figure  2  is  generalized  in  the  sense  that  it  shows  the  main 
components  of  the  process.  For  readers  who  are  unfamiliar  with  the  components  shown  in  Figure  2,  the 
following  descriptions  will  be  helpful: 
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GESAR:  An  operating  system  fw  the  TI4100.  GESAR  was  developed  by  the  Naval  Surface  Warfare  Center 
(NSWQ  and  is  loaded  into  the  TI4100  memory  when  in  use.  It  was  designed  to  collect  GPS  data 
with  geodetic  applications  in  mind.  Details  regarding  its  operation  and  outputs  can  be  found  in 
Darnell  [1987], 

BEPP/CQRE:  An  operating  system  for  the  TI4100  developed  by  the  Applied  Research  Laboratories  of  the 
University  of  Texas  at  Austin  (ARL/UT)  for  the  environment  where  the  receiver  is  driven  by  a 
personal  computer.  Such  an  environment  exists  at  the  five  DMA  GPS  Hxed  tracking  sites  which 
support  GPS  precise  orbit  determinadon.  The  BEPP  (Basic  External  Processor  Program)  resides  in 
the  laptop  PC  while  CORE  resides  in  the  TI4100. 

FICA:  An  acronym  for  Eixed,  Integer,  Character.  ASCII  formatted  data  files.  This  data  structure  was 
developed  by  ARLAJT  for  a  variety  of  receiver  types.  A  series  of  specific  block  idendfiers  has 
been  assigned  for  specific  uses.  For  example,  FICA  Block  6  contains  the  pseudorange  and  carrier 
beat  phase  data  from  the  four  trackers  of  the  714100  when  GESAR  is  used.  FICA  Block  SS 
contains  this  informadon  (with  addiuonal  quandties)  when  BEPP/CORE  is  used.  This  file 
structure  preserves  the  original  data  recorded  by  the  receiver  and  is  thus  suitable  for  archiving 
purposes.  It  is  also  suitable,  in  many  cases,  for  data  exchange  with  other  organizadons.  A 
non-ASCII  (machine  dependent  binary)  structure  (FIQ  is  also  available  if  desired.  A  complete 
descripdon  of  the  FICA  structure  can  be  found  in  ARL/UT  [1987]. 

EPHEMERIDES  AND  CLOCKS:  The  ’post-fit’  (precise)  GPS  ephemerides  and  satellite  clock  state 
esdmates  are  obtained  from  DMAHTC/GG  as  output  from  the  weekly  (8-day)  orbit  determinations 
using  the  OMNIS  program.  The  sateUite  posidon  and  velocity  vectm^  are  given  in  the  WGS  84 
reference  frame  at  15  minute  intervals  while  the  clock  state  estitnates  are  given  at  hourly  intervals. 
A  complete  mathemadcal  formuladon  of  OMNIS  is  found  in  Swift  [1987],  Note  that  the 
broadcast  ephemerides  and  clock  state  prediedons  are  available  in  the  RCA  (or  FIQ  file.  An 
assessment  of  the  weekly  precise  ephemerides  and  clock  states  has  been  presented  by  Gouldman  et. 
al.  [1989], 

STARPREP:  The  GEOSTAR  PREPROCESSOR  described  in  Section  6. 

POINT  FILE:  A  binary  file  produced  by  STARPREP  which  contains  identifying  information,  the  corrected 
tracking  data,  individual  correction  values,  and  interpolated  GPS  ephemerides  and  clock  states. 

GASP:  The  Qeodedc  Absolute  Sequential  Eosidoning  program  described  in  Section  7. 

Users  of  the  algorithms  and  programs  described  here  will  find  it  highly  desirable  to  perform  as 

many  steps  as  possible  on  the  same  computo'  equipment  Experience  has  shown  that  transfer  of  data  files 

from  one  machine  to  another  causes  the  most  delay  in  the  processing  flow. 
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2.  MISSION  PLANNING  AND  OBSERVATION 
REQUIREMENTS 

Since  the  TI4100  receiver  is  limited  to  collecting  data  from  4  satellites  at  a  time,  a  range  of 
possible  tracking  geometries  will  be  available.  Hie  user  must  design  a  data  collection  scheme  which  will 
maximize  the  geodetic  utility  of  the  data  set  Primarily,  this  involves  selecting  the  appropriate  hours  for 
tracking  and  the  specific  satellite  combinations  that  are  to  be  tracked.  This  process  is  referred  to  as 
selecting  scenarios.  Ideally,  a  geodetic  user  would  like  to  encounter  a  geometric  situation  where  GPS 
satellites  are  present  in  all  four  quadrants  of  his  local  hcxizon.  Moreover,  the  satellites  should  be  observed 
through  a  range  of  elevation  angles.  Sky-plots  such  as  those  shown  in  Figure  3  are  very  helpful  when 
performing  mission  planning.  They  allow  an  easy  visualization  of  azimuth  and  elevation  angle.  The 
tracking  geometry  during  any  particular  data  collection  period  will  primarily  be  a  function  of  geographic 
location.  The  GPS  orbital  characteristics  cause  the  current  (Block  I)  constellation  to  be  available,  at  any 
given  location,  approximately  4  minutes  earlier  each  day.  This  rule  of  thumb  can  be  used  to  avoid 
generating  mission  planning  plots  for  each  day  of  a  tracking  campaign.  The  distribution  of  the  Block  I 
satellites  purposely  optimizes  the  tracking  geometry  for  stations  located  in  the  southwestern  US.  Although 
plots  such  as  those  shown  in  Figure  3  are  helpful  during  the  planning  phases  of  a  survey,  one  should  be 
aware  of  the  limitations  of  such  plots.  For  example,  they  do  not  reflect  the  local  environment  where  the 
receiver  is  to  be  placed  and  thus  cannot  account  for  local  obstructions  such  as  mountains,  buildings,  heavy 
forestation,  power  lines,  towers,  fences,  etc.  The  Held  personnel  may  find  that  local  conditions  such  as 
these  severely  affect  the  planned  data  collection  scenario.  The  field  personnel  should  attempt,  if  possible,  to 
remove  all  obstructions  such  that  the  receiver  can  successfully  track  (with  minimal  multipath)  down  to  ten 
degrees  elevation  angle.  If  obstructions  caimot  be  cleared,  construction  of  a  tower  may  be  necessary  to  raise 
the  antenna  above  the  obstructions. 

Note  that  the  elevation  angles  shown  in  Figure  3  (and  similar  plots)  are  measured  from  the  tangent 
to  the  ellipsoid  (or  sphere)  and  thus  do  not  account  for  the  geoid  or  the  local  terrain. 

Experience  has  shown  that  the  quantity  known  as  Epsition  Qjlution  ofEfecision  (P-DOP)  can  be 
helpful  in  evaluating  the  geometric  situation.  Hiis  quantity,  along  with  a  subset  of  similar  quantities 
(G-DOP,H-DOP,  etc.)  is  available  from  some  survey  planning  packages.  For  point  positioning  as  well  as 
relative  positioning,  a  falling  P-DOP  is  highly  desirable.  The  tracking  schedule  should  always  be 
organized  around  the  times  of  falling  P-DOP. 

Since  the  TI4100  tracks  only  four  satellites  at  a  time,  and  more  than  four  satellites  may  be  in  the 
user's  sky  during  a  backing  session,  a  change  in  scenario  during  a  tracking  session  is  possible.  These 
scenario  changes  are  desirable  for  point  positioning  since  they  provide  an  expanded  range  of  geometries 
during  the  tracking  session.  Experience  gained  to  date  indicates  that  a  geodedc-quality  point  position  (a 
position  with  a  standard  deviation  on  each  component  of  less  than  one  meter)  can  easily  be  obtained  with 
the  GASP  algorithm  when  the  data  set  spans  a  period  longer  than  about  4  hours.  If,  however,  the 
geometric  situation  is  very  favorable,  a  geodetic-quality  point  position  can  be  achieved  with  GASP  from  a 
smaller  amount  of  data.  The  GASP  output  contains  a  number  of  quality  control  measures  which  the  user 
can  examine  to  evaluate  the  precision  of  the  estimated  position.  The  minimum  observation  requirements 
for  geodetic  point  positioning  are  summarized  in  Table  1. 

Since  the  differencing  algorithm  employed  in  GASP  removes  most  of  the  effects  of  the  receiver's 
clock,  there  is  no  requirement  to  utilize  an  external  frequency  standard  with  the  TI4100  receiver. 
Experiments  have  been  performed  which  compared  results  obtained  from  receivers  which  operated  on  quaru 
crystal  and  rubidium  oscillators.  No  difference  in  the  quality  of  positioning  results  was  detected.  This  issue 
will  have  to  be  reexamined  when  data  is  collected  in  the  AntiSpoofing  /  Selective  Availability  (AS/SA) 
environment.  All  developmental  testing  and  evaluation  of  the  GASP  and  STARPREP  algorithms  has  been 
done  with  the  data  coUection  interval  set  at  30  seconds. 

The  a  priori  (starting)  station  coordinates  which  STARPREP  and  GASP  use  are  usually  assumed 
to  have  standard  deviations  of  1/2  kilometer  in  each  component  In  general,  the  starting  coordinates  will  be 
more  precise  than  this.  These  coordinates  are  contained  in  the  FICA  file  (Block  101  with  GESAR  or  Block 
124  with  BEPP/CORE).  They  are  input  by  the  TI4100  operator  at  the  teginning  of  the  tracking  session. 
They  can  be  updated  or  edited  in  the  FICA  file  if  necessary. 
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Table  1.  Minimum  Observation  Requirements  for  Geodetic  Point  Positioning  with  the  TI41(X) 


Criteria  \  Desired  Precision 


Standard  Deviation  (1  Sigma  in  each  position  component) 
0.75  meters  1.5  meters 


Minimum  observation  period 

Data  collection  interval 

Minimum  number  of  S  Vs  to  be 
tracked  simultaneously 

Minimum  number  of  scenarios 
of  4  satellites 

Minimum  number  of  S  V  orbit 
planes 

PDOP 


Satellite  geometry 


TI4100  operating  system 

Frequency  standard 

Minimum  crystal  oscillator 
warm-up  time 

Meteorological  observations 


Antenna  heights  and  local 
survey  ties 


240  minutes  90  minutes 

30  seconds  30  seconds 


4  4 


2  1 


PDOP  for  each  scenario  is  to  fall  as  much  as  possible 
over  the  period  of  observation.  For  the  last  40  minutes, 
the  PDOP  should  be  less  than  10  and  no  greater  than  5  at 
the  end  of  the  tracking  session. 

The  SVs  should  be  distributed  in  all  four  quadrants  of  the 
observer's  sky.  At  least  one  SV  should  appear  within 
of  the  prime  vertical,  anywhere  within  the  elevation  angle 
range  of  10°  to  50°.  All  observations  are  to  be  above  10° 
elevation  angle. 

GESARorBEPP/CORE 

Internal  crystal  or  external  atomic 

1  hour  1  hour 


Recorded  at  commencement,  each  hour,  and  on  completion 
of  tracking  session.  Temperatures  are  to  be  accurate  to  tl°C 
barometric  pressure  to  ±1  mb  and  relative  humidity  to  ± 
5%. 

Height  from  the  monument  to  the  base  of  the  antenna  is  to 
determined  to  ±lcm.  All  local  survey  ties  should  also  be 
accurate  to  ±1  cm.  An  annotated  sketch  of  the  local  survey 
is  essential. 


3.  FIELD  VALIDATION 

Before  leaving  a  remote  tracking  site,  a  TI4100  user  should  know  that  the  data  collected  will  be 
geodetically  useful.  If  the  mission  planning  stage  is  performed  correctly,  the  tracking  session  followed  the 
mission  plan,  and  the  equipment  performed  normally,  there  would  be  no  need  for  any  in-field  validation. 
Of  course,  a  number  of  unexpected  events  can  occur  before  or  during  a  tracking  session  which  necessitate  a 
change  in  plans.  Such  events  could  include  a  satellite  becoming  'unhealthy',  a  delay  in  getting  personnel 
to  a  remote  tracking  site,  malfunctioning  equipment,  severe  environmental  conditions,  etc. 

Most  forms  of  in-field  validation  require  computer  resources  and  a  power  supply.  If  the  TI4100 
data  is  recorded  on  a  cassette  tape,  a  tape  reader  may  also  be  necessary.  If  resources  such  as  these  are 
available  in  the  field,  a  number  of  software  packages  can  be  utilized  to  verify  that  the  data  will  be 
geodetically  usefiil.  These  packages  include  but  are  not  limited  to: 


PROGRAM  NAME 

SOURCE 

a.GS2nCandnCnCA 

ARL:UT 

b.PRTN 

DMASC/EG 

c.  DSUMI'Y 

ARL:UT 

BASIS  FOR  VALIDATION 

Detects  certain  errors 
Provides  number  of  data  blocks  recorded 
ASCn  (FICA)  file  can  be  examined  directly 
Summarizes  all  data  collected 
Predicts  data  quality  using  GDOP,  PDOP 
Detects  anomalies  and  receiver  entas 
Summarizes  data  blocks,  satellites, 
observation  periods 


Each  of  the  above  computer  programs  is  operational  on  an  IBM  PC/XT  (or  compatible 
equipment).  In  the  event  that  computer  resources  are  not  available  in  the  field,  the  DSUMRY  program  can 
be  loaded  on  the  TI4100  memory  and  used  to  verify  the  data  file's  contents. 

In  all  cases,  the  field  personnel  must  record  accurate  field  notes  regarding  the  local  conditions  and 
report  any  anomalous  incidents  such  as  thunderstorms,  a  weather  front  passing  through,  electronic 
malfimctions,  power  failures,  etc.  The  following  information  is  essential  in  the  post-processing  stage: 

V  Receiver  operating  system  system  and  version  number  (ex:  GESAR  vl.5,  BEPP/CORE  v2.0) 

V  Type  of  fr^uency  standard  used  fintoTial  crystal;  external  rubidium,  cesium,  maser) 

V  Weather  observations;  Time,  Temperature  (®Q,  Pressure  (mb),  Relative  humidity(%) 

V  Accurate  antenna  height  measurements  (fiom  monument  to  base  of  antenna) 

V  Accurate  local  survey  tie  data,  if  applicable,  including  detailed  sketch 

V  The  TI41(X)  navigation  solution  and  any  other  reliable  a  priori  position  estimate 

Even  when  information  such  as  this  is  recorded  in  the  data  set  itself,  a  set  of  accurate  field  notes  is 
to  be  maintained  which  clearly  record  this  important  information.  Standard  field  log  forms  have  been 
developed  for  this  purpose.  A  complete  set  of  field  log  forms  is  given  in  Appendix  A. 
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4.  TRANSLATION  OF  GESAR  OR  BEPP/CORE  TO  THE  FICA 
STRUCTURE 

4.1  Overview 

If  the  tracking  data  collected  with  the  GESAR  or  BEPP/CORE  operating  systems  is  to  be  used  for 
GASP  point  positioning,  it  must  be  translated  into  the  FICA  file  structure.  This  is  the  structure  required 
by  the  STARPREP  preprocessor.  If  the  original  data  is  a  GESAR  file  on  cassette  tape,  a  minimum  of 
three  utility  programs  must  be  executed  in  order  to  accomplish  this  translation.  These  programs  are  named 
MFERD,  GS2FIC  (GESAR  to  HC  or  BEPP/CORE  to  HC)  and  FICHCA  (HC  to  HCA).  These 
uanslation  programs  have  been  developed  by  ARL/UT.  If  the  original  data  file  was  generated  with 
BEPP/CORE,  the  file  will  be  on  a  diskette.  In  this  case,  only  the  GS2FIC  and  FICFICA  utilities  are  used. 
Versions  of  these  programs  for  a  personal  computer  are  available.  As  mentioned  above  (section  3).  this 
translation  process  can  serve  as  a  type  of  field  validation  if  it  is  performed  before  the  data  set  is  sent  for 
post-processing.  Upon  execution  of  the  GS2FIC  program,  a  tabulated  summary  of  HC  block  types  is 
presented.  A  detailed  set  of  standard  operating  procedures  for  this  translation  and  electronic  transfer  of  the 
data  files  is  given  below. 


4.2  Data  Translation:  Operating  Procedures 


This  section  describes  a  routine  procedure  for  transforming  TI4100  GESAR  or  BEPP/CORE 
tracking  data  from  cassette  tape  (GESAR)  or  diskette  (BEPP/CORE)  to  9-track  magnetic  tape  for  transport 
to  a  computer  where  the  point  positioning  program  (GASP)  and  its  preprocessor  (STARPREP)  reside. 
These  procedures  were  develop^  from  Ison  [1988]  and  Yambot  [1988],  and  are  specific  to  the  ^^-DOS 
based  ^nith  PC  and  the  DEC  MicroVAX  computer  equipment  at  the  DMAHTC/GGSC  facility. 

The  FICA  file  structure  requires  about  S  times  more  storage  space  than  the  original  binary  data 
(GESAR  or  BEPP/CORE).  The  FICA  structure  facilitates  portability  since  these  ASCII  files  are  machine 
independent  The  STARPREP  pieprocesstv  was  designed  to  read  these  FICA  files.  If  a  continuous  data  set 
is  stored  as  a  series  of  files  on  the  Zenith  PC,  they  should  be  concatenated  in  chronological  order  before 
conversion  to  FIC  and  FICA  structures.  With  the  Block  I  GPS  constellation,  this  usually  only  occurs  with 
BEPP/CORE  data  from  the  DMA  fixed  sites. 

Note:  To  avoid  confusion  on  the  instructions  of  the  remainder  of  this  section  (4.2),  the 
responses  by  the  operator  will  be  shown  in  bold&ce  type. 

4J.1  Reading  GESAR  Data  From  the  Cassette  Tape  to  the  21enith  PC 


a.  Connect  the  interface  cable  from  the  MEJ^fTEC  5450XL  reader  to  the  back  panel  of  the  Zenith  PC.  The 
port  at  the  Zenith  PC  is  the  first  one  from  the  left  viewed  from  behind  (pcvts  and  connectors  should  be 
properly  labeled). 

b.  On  the  surge  suppressor  box  (or  the  power  protector  box)  turn  the  power  switch  on. 

c.  On  surge  suppressor  box  only,  press  the  Reset  button. 

d.  On  the  MEMTEC  cassette  upe  reader,  turn  power  on. 

e.  Load  data  cassette  tape  into  MEMTEC  tape  reader.  Make  sure  the  cassette  tape  is  at  the  beginning. 

If  not,  press  the  "RWND"  button  on  the  MEMTEC  reader  to  rewind  it  to  the  beginning. 

f.  On  the  Zenith  PC,  access  the  D  disk  drive  by  typing  in:  D:  <Return> 

g.  To  change  the  directory  to  the  FICA  directory,  type  in:  CD  FICA  <Return> 
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h.  Execute  the  program  that  reads  the  cassette  tape  by  typing:  MFERD  <Retuin> 

i.  The  program  prompts  for  the  output  file  name.  Type  in  a  file  name  of  twelve  charactws  in  the 
following  way; 

Character  I,  "G"  indicates  GESAR  data,  "B"  indicates  BEPP/CORE  data 

Character  2,  ”B”  indicates  binary  image  file 

Character  3-6,  last  four  digits  of  station  ID 

Character  7,  cassette  number  or  session  letter 

Character  8,  last  digit  of  the  year 

Character  9, period  required  by  MS-DOS 

Character  10-12,  Day  number  of  year,  when  data  was  collected 

For  example,  the  file  GB5271A9.148,  will  represent  GESAR  binary  data  from  station  85271,  cassette 
session  A,  of  year  1989  on  day  148. 

(Enter  input  file  name)  <Retum> 

j.  Now  the  program  asks  for  the  type  of  dump.  Zero  (0)  for  binary  and  one  (1)  for  hexadecimal.  Select 

binary  by  typing  in:  0  <Return> 

k.  Press  "Load  Point"  button  on  the  cassette  reader. 

l.  Press  the  "Return"  key. 

m.  Reading  process  starts  when  the  "Busy",  "Data",  and  "Ready"  indicators  flash  and  "BIN”  button 
flashes  on  the  MEMTEC  reader.  When  reading  is  completed  the  prompt:  "Stop  -  Program  terminated" 
appears.  After  this  message  appears,  press  the  "RWND"  button  to  rewind  the  cassette  tape. 

n.  Remove  the  cassette  tape. 

o.  Check  the  new  file  name  and  its  number  of  bytes  on  the  disk  directtvy  by  typing:  DIR  <Return> 

p.  Log  the  file  name  and  the  number  of  bytes  on  a  processing  log  sheet 


4  J,2  Creating  the  FIC  file 

The  GESAR  or  BEPP/CORE  binary  files  may  be  translated  to  FIC  and  FICA  on  either  the  Zenith 
PC  (w  the  DEC  MicroVAX.  If  the  latter  is  used,  the  data  is  transferred  from  the  PC  to  the  MicroVAX 
using  the  procedures  in  section  4.2.4.  The  FIC  file  is  then  created  following  steps  c  through  i  in  this 
section.  The  MicroVAX  command  procedures  are  located  in  the  directory:  [GSGC.HTLRM. 
EXPORT.COMMAND] 

a.  If  the  current  disk  is  not  the  D  drive,  type  in:  D:  <Return> 

b.  Access  the  FTCA  directory  by  typing:  CD  FICA  <Retum> 

c.  To  execute  the  necessary  program,  type  in:  GS2FIC  <Return>  (@GS2FIC  <Retuni>  on  MicroVAX) 

d.  The  program  prompts  for  the  "INPUT  FILE  NAME".  This  is  the  file  that  was  read  from  a  cassette  tape 
(GESAR)  or  diskette  (BEPP/CORE)  onto  the  Zenith  PC.  (Enter  input  filename)  <Return> 

(on  the  MicroVAX,  the  file  name  is:  [GSGCJ>EPOTJFICA]filename  <Retum> ) 

e.  ■rae  program  now  prompts  you  for  the  "OUTPUT  FILE  NAME”.  A  maximum  of  twelve  characters  is 

allowed.  Use  the  naming  convention  specified  in  Section  4.2.1  except  on  Character  2  where  "B"  will 
become  "F”  to  define  a  FIC  output  file.  (Enter  output  filename)  <Retum> 

(on  the  MicroVAX,  the  file  name  is:  [GSGCJIEPOT  JlCAlfilename  <Retum> ) 

f.  Next,  the  program  prompts  for  "ORBIT  FILE  NAME".  Type  in;  NUL  <Renmi> 

g.  Then,  the  program  asks  for  header  information  with  the  message  "ENTER  TITLE  UP  TO  37 
CHARACJTERS".  This  header  will  become  the  first  record  in  the  FIC  and  FICA  file.  This  step  is 
critical  since  this  record  is  the  only  means  of  identifying  the  file's  contents  once  it  is  put  on  a  9-Track 
tape.  Enter  the  information  in  the  following  way; 

Character  1  -9,  Project  Name; 

Character  10-19,  Station  Number  or  Station  Name; 

Character  20-29,  Year, 

Character  30-37,  Day  of  year,  when  data  collected  started; 


It  is  left  justified.  Leave  blanks  in  remaining  spaces.  Hoe  is  an  example: 

GREENLAND81399  1989  176 

Press  the  "Return"  key. 

h.  While  the  program  is  running,  messages  describing  skipped  data  will  be  displayed.  These  error 
messages  are,  in  most  cases,  normal  and  are  the  result  of  the  method  by  which  data  is  recorded.  When 
finished,  several  messages  appear  on  screen .  One  of  them  states  "GS  2  FIC  CONVERTED". 

i.  Log  the  created  file  name  and  number  of  bytes  on  the  file  management  sheet  (a  sample  of  a  file 
management  sheet  is  given  in  Appendix  A). 

MdDiriSr  Files  in  the  EI£  format  cannot  ts  transferred  from  one  machine  to  another  since  this  structure  is 
machine-dependent  binary. 

4.2  J  Creating  the  FTCA  file 

The  computer  which  created  the  FIC  file  must  be  used  to  create  the  FIC  A  file.  If  the  MicroVAX 
was  used  to  create  the  FIC  file,  follow  steps  c  through  f  in  this  section  to  create  the  FICA  file.  The 
MicroVAX  command  procedures  are  located  in  the  directory:  [GSGC.IirLRM£XPORT.COMMAND] 

a.  If  the  current  disk  is  not  the  D  disk  drive.  Type  in;  D:  <Retum> 

b.  Change  into  the  FICA  directory  by  typing:  CD  FICA  <Retum> 

c.  Execute  the  program  by  typing:  FICRCA  <Retutn>  (@FICFICA  <Rctiim>on  the  MicroVAX) 

d.  The  program  prompts  you  for  the  input  file  name.  This  is  the  FIC  file  that  was  created  in  an  earlier 
process  (see  Section  4.2.2).  (Enter  input  filename)  <Rctum>  (on  the  MicroVAX,  the  file  name 
is:  tGSGC.DEPOT.nCAirilename  <Retum> ) 

e.  The  program  now  prompts  for  the  output  file  name.  Again,  a  maximum  number  of  twelve  characters  is 
allowed.  The  input  and  output  file  names  should  be  identical  except  on  Character  2  where  "P  (for  £IC) 
will  become  "A"  (for  FIC^.  The  "A"  represents  ASCII.  (Enter  output  filename)  <Retum> 

Q  n  the  MicroVAX,  the  file  name  is:  [GSGCJIEPOT JICAlfilename  <Retum>  ) 

f.  When  finished:  "FIC  TO  FICA:  TRANSLATED  ##  RECORDS"  will  appear  on  the  screen.  Here, 

"##"  represents  an  integer  number  that  is  displayed  and  it  is  the  total  number  of  translated  records. 

Log  the  file  name  and  number  of  records  on  the  file  management  sheet 

4.2.4.  Transferring  GESAR,  BEPP/CORE  or  FICA  Files  from  the  Zenith  to  the 
MicroVAX 

a.  Connect  the  gray  interface  cable  from  the  MicroVAX  computer  to  the  first  port  from  the  left  as  viewed 
from  behind  on  the  back  panel  of  the  Zenith  PC. 

b.  Set  the  directory  where  the  files  to  be  transferred  are  located.  This  will  usually  be  the  FICA  directory  on 
the  D  disk.  The  transfer  program  is  KERMIT.  To  execute  this  program,  type  in  this  name  and  a 
pre-defined  directory  path  will  take  you  directly  into  KERMIT.  With  any  prompt,  type  in: 

KERMIT  <Returo> 

c.  The  KERMIT  prompt  will  appear.  Connect  to  the  MicroVAX  computer  by  typing: 

Kefmit-MS>CONNECT  <Rctum> 

d.  Hit  the  Return  key  two  or  three  times.  The  screen  will  clear. 

e.  The  MicroVAX  computer  will  ask  you  for  your  username  and  password.  Type  these  in: 

Usentame:  GSGC  <Return> 

Password:  (Enter  your  password  here)  <Retum> 

f.  After  access  is  gained  to  the  MicroVAX  system,  the  $  prompt  will  appear.  Type  and  enter  the  directory 

name  where  the  file  is  to  be  sent  Enter  $  FICA  (dir^tory-name)  <Retum> 


g.  Next,  enter  $  KERMIT  <Retum> 

h.  If  file  to  be  transferred  is  ASCII  (an  FICA  file)  type,  slop  to  step  i.  If  the  file  is  binary  (a  GESAR  or 
BEPP/CORE  file),  then  type  and  enter  Kermit-32>SET  FILE  TYPE  BINARY  <Retum> 

MtOTlS:  FIC  files  cannot  be  transferred;  only  GESAR,  BEPP/CORE  or  FICA  can. 

i.  Type  and  enter  Keiniit-32>SERVER  <Retum> 

j.  A  long  message  will  appear  on  bottom  of  the  screen.  Depress  the  CTRL  and  the  right  bracket  I  keys 

simultaneously  and  then  release  them.  Then  press  and  release  the  C  key: 

k.  Send  the  file  of  your  choice  from  the  Zenith  to  the  MicroVAX  by  typing: 

Kermit-MS>SEND  D:\FICA\filename  <Retum> 

Here,  "filename"  is  the  name  of  the  file  being  sent  Immediately  the  screen  is  cleared,  and  transfer 
information  is  displayed.  For  example,  if  the  file  GB5271A9.148  is  being  transferred,  the  top  of  the  screen 
will  appear  as  follows: 

File  name:  GB5271A9.148 
KBytes  transferred:  3 
Percent  Transferred;  25% 

Sending:  In  progress 

Some  of  the  numbers  will  be  changing  fast  and  continuously.  When 

"Percent  transferred"  reads  100%  and  "Sending"  reads  Completed,  the  transferring  is  finished. 

If  another  file<s)  is  to  be  transferred,  repeat  this  step  (}c)  with  the  new  file  name(s). 

l.  Then  type:  Kennit-MS>FINISH  <Retum> 

m.  Now,  just  type:  Kennit-MS>C  <Retum> 

n.  Press  the  "Return"  key  to  connect  back  to  the  MicroVAX  computer.  The  ’$’  prompt  will  appear. 

o.  To  verify  a  successful  transfer,  check  the  file  name(s)  in  the  directory.  Type  and  enter 

$  DIRECTORY  <Retiim> 

p.  If  a  FICA  file  was  transferred,  verify  the  first  few  Blocks  in  the  transferred  FICA  file  by  viewing  it  with 

the  editor  $edit  [GSGC.DEPOT.FICA]filenaine 

Exit  the  MicroVAX  editor  by  typing  SHIFT  and  PFl  keys  simultaneously,  then  press 

Q 

q.  If  transferring  is  complete,  just  repeat  steps  g,  i,  j,  and  exit  from  Kermit  by  entering: 

Kennit-MS>BYE  <Retum> 

This  will  return  control  to  the  Zenith  PC  in  the  directory:  D:NFICA 


4,2 J  Saving  a  Blocked  FICA  File  Onto  9'Track  Tape  From  the  MicroVAX 

a.  This  procedure  can  be  performed  either  from  the  Zenith  (as  a  MicroVAX  terminal)  or  the  MicroVAX 
computer.  If  the  Zenith  is  used,  follow  steps  a  through  e  of  Section  4.2.4  to  connect  to  the  MicroVAX. 
If  using  the  MicroVAX  directly,  follow  from  step  e  in  that  same  section. 

b.  Set  the  default  directory  as  follows:  $SET  DEF  $DISK2;[GSGC.DEPOTJFICAl 

c.  If  the  file  name  needs  to  be  changed,  type  and  enter  the  following: 

S  RENAME  (Old  filename)  (New  niename)<Retum> 

d.  To  get  the  total  number  of  records  of  the  file,  call  the  line  counter  program  by  typing: 

$REC  (niename)<Retum> 

e.  Repeat  step  d  until  all  the  number  of  records  for  all  of  files  are  logged. 

f.  The  files  can  now  be  written  on  9-track  tape.  Load  a  blank  tape  with  write  ring  on,  into  the  Digital 
TS05  Tape  Drive  (press  the  "Load  and  "On-Line"  soft  padded  buttons). 


g.  Execute  the  program  to  save  and  block  the  flies  by  typing  this: 

S  @[GSGCJ>EPOT]100BLOCK.COM  <Return> 

This  program  will  write  ASCII  records  (80  characters  long)  in  100  record  tape  blocks  until  the  entire  file 
has  been  written  U3  tape.  An  end  of  file  mark  will  be  placed  at  the  end  of  the  file. 

h.  The  program  will  prompt  for  the  file  name  to  be  saved  on  tape.  Enter  the  file  name  you  want  Note  that 
the  default  directory  has  already  been  set  in  step  b  so  the  entire  path  name  need  not  be  specified. 

$  (Enter  filename)  <jletum> 

i.  At  the  message  "MOUNT  THE  TAPE&  PRESS  <RETURN>  TO  CONTINUE:",  press  the  Return 
key. 

j.  At  the  message  "NUMBER  OF  FILES  TO  SKIP?”  enter  the  aiqvopriate  number  of  files  to  skip  over 
which  may  have  already  been  written  on  the  tape.  If  the  tape  was  bl^  and  this  is  the  first  file  teing 
written  this  question  will  not  be  asked.  If  you  are  writing  over  a  previously  used  tape,  enter  0.  for  the 
second  file  enter  1 ,  for  third  enter  2,  etc.  Press  the  "Return"  key  after  your  answer.  The  file  will  then 
be  written  to  the  tape. 

k.  When  finished,  the  screen  will  display  the  number  of  input  records  and  the  number  of  output  records  (or 
the  number  of  blocks).  Log  these  numbers.  Compare  number  of  input  records  here  with  the  one 
obtained  from  the  line  counter  program  for  the  same  file.  Make  sure  they  are  the  same.  If  they  are  not, 
check  the  tape,  the  tape  drive,  or  get  help  from  the  MicroVAX  room  personnel. 

l.  Upon  completion,  the  following  message  is  displayed  on  the  screen: 

DO  YOU  WANT  TO  COPY  ANOTHER  FEi?  <YES  OR  NO>. 

If  other  files  are  to  be  written,  answer  YES  and  repeat  steps  h  through  I.  If  no  further  files  are  to  be 
copied  answer  NO.  The  tape  will  be  rewound  and  unload  automatically.  Press  the  "Return”  key 
after  all  answers. 

m.  Remove  the  tape  from  the  tape  drive  when  finished  with  all  files  and  log  oE  from  the  MicroVAX  by 

typing:  $  LO  <Retuni> 


4  J  FICA  Files  From  9-Track  Tapes 

Once  the  translation  is  successfully  performed,  the  FICA  file  may  have  to  be  transported  to  another 
computer.  Typically,  the  9'track  tape  generated  for  this  transfer  will  also  be  used  for  data  archiving.  The 
recommended  tape  density  is  16()0BPI  since  this  capability  is  common  to  most  equipment.  If  the  data  is 
exchanged  with  another  organization  at  a  later  date,  regenerating  the  tape  will  not  be  necessary  (a  simple 
ASCII  tape  copy  utility  can  be  used).  Experience  has  shown  that  30  to  40  FICA  files  (each  containing  6 
hours  of  30  second  data)  can  fit  on  one  2400ft  9-track  tape  when  1600BPI  is  used  and  the  tape  is  blocked  at 
100  records  per  tape  block.  Each  record  contains  80  ASCII  characters.  Simple  utility  programs  for  writing 
and  reading  blocked  ASCII  files  are  necessary  for  this  task.  Once  a  FICA  file  is  loaded  into  a  mass  storage 
file  from  the  9-track  tape,  it  can  be  treated  the  same  as  any  other  FICA  file. 

A  FICA  file  naming  convention  has  been  established  within  DMA.  The  structure  of  the  FICA  file 
name  is  as  follows: 


Example:  GEOSTAR*GA271A9148. 


Where  'GEOSTAR*'  is  the  Unisys  qualiHer  and  thus  would  not  be  necessary  on  non-Unisys 
equipment  The  remainder  of  the  file  name  is  the  critical  part 

Character  1  'G'  indicates  Gesar, 'B' indicates  BEPP/CORE 

indicates  ASCII  file 
indicates  last  3  digits  of  station  ID 
cassette  or  data  subset  sequencing  A,B,C.  etc. 
last  digit  of  year  data  was  collected 
day  number  of  year  that  data  was  collected 
use  for  edited  file  if  the  file  is  a  subset  of  the  original 
On  DEC  (VAX  or  MicroVAX)  equipment,  the  file  name  extension  '.DAT  should  be  used. 


1 

G' 

2 

'A' 

3-5 

771’ 

6 

'A' 

7 

'9' 

8-10 

11 

'148' 

13 


5.  SCREENING  AND  EDITING  FICA  FILES 


Once  a  FICA  file  is  loaded  onto  the  computer  where  it  will  be  processed,  it  can  be  examined  and 
manipulated  in  a  number  of  different  ways.  The  simplest  tool  one  can  use  is  any  appropriate  ASCII  file 
editor.  An  analyst  can  simply  browse  through  the  file  and  visually  verify  that  the  file  contains  the  properly 
structured,  necessary  information.  The  analyst  must  be  familiar  with  the  FICA  block  identifiers.  For 
example,  when  working  with  a  BEPP/CORE  -  created  file,  the  Block  lO’s  should  be  checked  to  verify  that 
the  weather  observations  in  the  data  file  match  those  recorded  on  the  field  sheet  If  no  Block  lO's  are  present 
(the  case  when  GESAR  is  used),  the  weather  data  from  Block  101  will  be  used  for  the  entire  data  set  If 
necessary,  the  editor  can  be  used  to  modify  the  ccMitents  of  the  file. 

Since  manually  screening  a  data  file  is  very  tedious  and  provides  only  a  limited  amount  of 
information,  a  comprehensive  screening  program  has  b^  written  at  DMA  to  allow  analysis  and  editing  of 
FICA  files  and  provide  a  summary  of  a  data  file's  contents.  This  program  has  been  designated  PRTN, 
which  is  a  contraction  for  PaRTitioN  a  data  file  (s).  This  utility  program  can  screen  and  partition  an 
individual  data  set  or  a  pair  of  simultaneous  data  sets  which  will  be  used  for  baseline  estimation  programs. 
A  user's  guide  for  PRTN  has  been  prepared  by  Ayrandjian  [1989].  The  point  positioning  user  may  desire  to 
split  up  a  data  set  before  processing,  if  for  example,  a  significant  part  of  a  d^  set  is  to  be  discarded  (if  it 
contains  only  2  satellites  for  example). 

PRTN  provides  a  comprehensive  summary  of  a  FICA  file's  contents  and  allows  selective  editing 
of  a  file(s)  by  choosing  options  from  a  menu.  One  essential  piece  of  information  obtained  from  PRTN  is 
a  list  of  satellites  which  are  present  in  the  data  set  being  analyzed.  The  user  will  always  need  to  know 
which  satellites  are  present  before  giving  the  data  to  STARPREP.  Moreover,  if  the  FICA  file  contains  less 
than  4  satellites  for  significant  periods  during  the  data  span,  these  periods  may  be  removed  from  the  data  set 
Whenever  PRTN  performs  any  editing  of  an  FICA  file,  it  maintains  the  (»iginal  FICA  file  contents  and 
creates  a  new  FICA  file  which  will  contain  the  updated  contents.  The  label  TJPD'  is  added  to  the  new 
(edited)  FICA  file  name.  The  placement  of  this  identifier  dqrends  on  which  computer  equipment  PRTN  is 
being  used  on.  If  PRTN  is  used  to  modify  the  contents  of  a  file,  the  updated  FICA  file  is  to  be  used  from 
that  point  on. 

Independent  of  PRTN,  another  utility  program,  named  EDITFICA,  may  also  be  used  to  divide  a 
FICA  file  into  subsets.  This  relatively  simple  program  prompts  the  operate  for  the  start  and  end  tiroes  of 
the  data  subset  being  created.  This  program  is  especially  useful  when  operating  on  BEPP/CORE,  DMA 
fixed-site  data  sets  since  it  filters  out  the  numerous  weather  blocks  (Block  10s)  and  keeps  only  one  weather 
block  per  hour.  STARPREP  can  handle  up  to  five  weather  blodts  in  any  single  FICA  file. 

When  analyzing  the  contents  of  a  FICA  file  ftH-  positioning  use,  the  following  items  must  be 

checked: 

1.  A  valid  antenna  height  (From  Block  101  M’  Block  124).  This  antenna  height  must  be  the 
measured  distance  from  the  monument  to  the  base  of  the  TI4100  antenna.  GASP  applies  the  calibrated 
distance  from  the  base  of  the  antenna  to  the  mean  electrical  center  of  the  TI  antetma.  The  value  in  the  FICA 
file  should  agree  with  the  value  written  on  the  field  sheets.  Note  that  this  antenna  height  may  change  from 
day  to  day  if  the  receiver  occupied  a  site  for  repeated  days. 

2.  Valid  weather  data  (From  Block  101  Of  Block  10)  If  default  values  are  used,  a  degradation 
in  positioning  results  will  occur  which  is  difficult  to  detect  by  die  user.  If  default  values  are  found,  or  any 
error  is  detected  in  the  values  in  the  FICA  file,  the  user  should  correct  these  errors  befme  continuing  with 
processing.  The  weather  values  in  the  FICA  file  should  be  checked  against  the  values  written  on  the  field 
sheets.  The  barometric  pressure  values  can  be  checked  for  gross  errors  by  comparing  to  a  crude  value  for 
atmospheric  pressure  which  can  be  estimated  from  the  station's  elevation  above  mean  sea  level.  Charts  are 
available  for  this  conversion.  If  a  gross  error  is  detected,  it  may  be  due  to  the  use  of  a  barometer  that  has 
not  been  calibrated.  If  this  occurs,  the  user  of  the  data  should  not  resort  to  a  default  atmospheric  pressure. 
Instead,  the  user  should  use  a  value  implied  by  the  station's  elevation  above  mean  sea  level. 

When  the  GESAR  operating  system  has  been  used,  one  set  of  weather  values  will  be  found  in 
Block  101  of  the  FICA  file.  When  BEPP/CORE  has  been  used  to  record  data,  a  series  of  weather  Blocks 
(Block  10)  will  be  present  in  the  data. 
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3.  The  year  and  the  dav  number  (Block  101  or  Block  124)  that  the  data  file  was  recorded 
should  be  noted  so  the  user  can  assign  the  appropriate  precise  ephemeris  files  when  the  preprocessor  is  run. 
If  the  FICA  file  contains  data  after  12:00  UT  on  a  Saturday  (S61600  seconds  into  the  GPS  week), 
'HALF- WEEK'  precise  ^hemeris  and  clock  files  must  be  used.  Tliese  HALF- WEEK'  precise  ephemeris  and 
clock  files  are  specifically  designed  to  allow  processing  of  data  which  was  collected  through  a 
week-crossover.  TTie  standard  precise  ephemeris  files  distributed  by  DMAHTC  consist  of  one  full  GPS  week 
(00:(X):(X)  Sunday  to  24:00:00  Saturday,  plus  Ifl  day  overlap  on  each  end).  The  "HALF-WEEK’  precise 
ephemeris  files  ^gin  at  00:(X):00  Wedhiesday  and  end  at  24:00:00  on  Tuesday.  There  is  no  smoothing 
performed  across  the  two  separate  GPS  weeks  when  the  HALF-WEEK  files  are  created. 

4.  The  user  should  take  note  of  any  week-crossovers  (Block  101  or  Block  124)  in  the  data  files. 
If  a  satellite  other  than  the  first  available  satellite  in  the  file  experiences  a  week  crossover  and  the  first 
available  satellite  data  span  falls  completely  into  the  later  of  the  two  weeks,  it  is  necessary  to  edit  the  file 
such  that  the  edited  file  tegins  in  the  later  of  the  two  wedcs.  The  day  number  in  the  Block  101  of  the  edited 
file  must  be  modified  to  represent  the  first  day  of  die  later  week.  If  the  first  satellite  tracked  experiences  a 
week  crossover,  there  is  no  need  to  edit  the  file  CHALF-WEEK'  precise  ephemeris  and  clock  files  must  be 
used). 

5.  The  a  oriori  station  coordinates  (Block  101  or  Block  124)  should  be  verified.  The 
nominal  accuracy  of  these  coordinates  is  assumed  to  be  half  of  a  kilometer  in  each  component  If  GASP 
rejects  a  large  percentage  of  data  on  the  basis  of  the  RMS  edit  or  if  the  residual  plots  produced  by  GASP  do 
not  appear  random  on  the  third  iteration,  there  is  probably  a  gross  error  in  these  a  priori  coordinates  in  the 
FICA  file.  If  improved  coordinates  are  not  available,  STARPREP  and  GASP  should  be  run  to  completion 
and  the  resulting  position  estimate  should  be  put  into  the  FICA  file  as  new  a  priori  coordinates.  Note  that 
STARPREP  must  be  run  again  (with  the  improved  a  priori  coordinates)  because  some  of  the  data 
corrections  arc  a  function  of  these  coordinates.  If  the  a  priwi  coordinates  ate  too  far  off,  some  of  the  data 
correction  applied  by  STARPREP  will  be  biased. 

6.  The  PRN  numbers  of  the  satellite  tracked  (Block  lOl  or  Block  124)  and  the  number 
of  satellites  in  the  data  set  should  be  noted.  Moreover,  any  data  gaps  or  periods  where  less  than  4 
satellites  have  been  tracked  should  be  noted.  The  length  of  the  data  span  should  also  be  noted. 

7.  The  following  FICA  Blocks  must  be  present  in  the  data  file  for  STARPREP  and  GASP  to 
produce  a  valid  point  position  estimate: 

G£SAB  DATA  BEPP/CORE  DATA 

Block  101  (Configuration)  Block  124  (Configuration) 

Block  9  (For  Broadcast  Ephemeris)  Block  9  (Fix  Broadcast  Ephemeris) 

Block  6  (Data  B  locks)  Block  55  (Data  Blocks) 

Block  10  (Weather) 

The  configuration  Block  (101  or  124)  must  appear  before  any  data  blocks. 

The  first  weather  block  (BEPP/CORE  only)  must  appear  before  any  data  blocks. 
Any  FICA  Blocks  not  shown  above  are  ignored  by  STARPREP  and  GASP. 

After  verifying  these  essential  pieces  of  information  discussed  here,  the  FICA  file  is  ready  for 
preprocessing  in  STARPREP.  If  the  RCA  file  was  edited  in  any  way,  the  new  (edited)  file  must  be 
checked  for  these  essential  elements  before  proceeding.  If  more  than  one  configuration  Block  (101  or  124) 
exists  in  a  file,  only  the  first  one  will  be  us^  by  STARPREP. 


One 

One  or  more 
Many 

One  or  more 
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6.  PREPROCESSING  (STARPREP) 

6.1  Overview 

The  GEOSTAR  PREPROCESSOR  program  is  designed  to  accomplish  a  number  of  essential, 
routine  iiincdons  which  are  necessary  to  exploit  the  data  for  the  maximum  precision  possible.  A  list  of  the 
standard  data  corrections  applied  by  STARPREP  and  their  functional  dependencies  is  shown  in  Table  2.  A 
detailed  mathematical  formulation  of  these  corrections  can  be  found  in  Malys  and  Ortiz  [op.cit].  The  main 
function  of  STARPREP  is  to  create  the  TOINT  File'  which  GASP  uses  to  perform  point  positioning.  The 
POINT  File  contains  the  corrected  LI  tracking  data,  the  interpolated  satellite  ephemerides  and  clock  states 
and  the  values  for  the  individual  data  correctitms  which  were  ^lied.  It  alM  contains  a  header  which 
identifies  the  data  set  Sample  STARPREP  runstreams  are  given  in  Section  6.2. 

TABLE  2.  Data  Corrections  and  Their  Functional  Dqtendencies 


Correction  /  Function  of 

Data 

Values 

Surface 

Weather 

Obs. 

A  priori 
Position 
Estimate 

Adopted 
Ephemeridej/ 
Clock  States 

1.  Epoch  of  Receipt  to 

Epoch  of  Transmission 

V 

2.  S^Uite  Clocks 

V 

3.  Receiver  Offsets* 

4.  Ionosphere 

V 

S.  Troposphere 

V 

V 

V 

6.  Generid  Relativity 

V 

7.  Earth  Rotation 

V 

>/ 

8.  Satellite  Antenna  Offsets 

V 

*  The  TI4100  frequency  offsets  of  -6000  Hz  and  4-7600  Hz  are  modeled  away  by  adding  6000(At)  cycles  and 
subtracting  7600(At)  cycles  to/from  the  L]  and  L2  recorded  carrier  beat  phase  respectively,  where  At 
represents  the  time  interval  (seconds)  between  the  data  value  epoch  and  the  current  reference  epoch  for  a 
particular  satellite. 


6.2  Sample  STARPREP  Runstreams 

Sample  runstreams  which  specify  the  selection  of  broadcast  or  precise  ephemerides/clocks  are  given 
in  Ingure  4.  The  following  7  data  corrections  are  to  be  ai^lied  for  all  geodetic  positioning  £q)plications: 

Range  Ionosphere  (IR)  Satellite  Clock  (SC) 

Doppler  Ionosphere  (ID)  General  Relativity  (GR) 

Ea^  Rotation  (ER)  Satellite  Antenna  Offset  (SA) 

Hopfield  Troposphere  (TH) 

These  corrections  are  selected  in  the  STARPREP  runstream  by  the  use  of  their  two-letter  codes 
shown  above.  The  'CMPUTCOR'  card  image  in  the  runstream  designates  which  corrections  are  to  be 
computed,  while  the  'APPLYCOR'  card  image  designates  which  corrections  are  to  be  applied  to  the 
pseudorange  and  converted,  carrier  beat  phase  data.  The  'PLOTCORR'  card  image  in  the  STARPREP 
runstream  is  used  to  designate  which  corrections  are  to  be  plotted. 
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The  cotiections  for  the  TI4100  receiver  frequency  offsets  are  automatically  applied  in  conjunction 
with  the  doppler  ionosphere  correction.  Further  details  regarding  the  frequency  offset  correction  are  given  in 
Malys  and  Jensen  [op.cit.],  and  in  the  note  at  the  bottom  of  Table  2. 

Occasionally,  there  may  be  reason  to  believe  that  a  particular  GPS  satellite  was  behaving 
erratically  or  experienced  some  anomaly  during  the  collection  of  a  givoi  data  set.  By  modifying  the 
STARPREP  runstream,  any  particular  satellite(s)  may  be  excluded  from  the  point  positioning  process. 
Satellites  which  operated  on  a  quartz  crystal  frequency  standard  during  the  data  span 
must  always  be  excluded.  To  accomplish  this  exclusion,  the  user  removes  all  lines  in  the  runstream 
which  refer  to  the  PRN  number  of  the  satellite(s)  to  exclude.  The  unit  numbers  for  the  file  assignments 
must  be  modified  such  that  they  are  continuous,  sequential,  in  ascending  order  and  commence  at  1 1  and  51 
(for  broadcast  preprocessing)  or  1 1, 21.  and  51  (for  precise  preprocessing).  The  excluded  sateilite(s)  must 
also  be  remov^  from  the  satellite  selection  list  represented  by  the  'SELECTTSV'  card  image.  The  sequence 
of  PRN  numbers  in  the  'SELECTS  V  card  image  must  match  the  sequences  of  PRN  numbers  in  the  file 
assignments.  Examples  are  given  in  Figure  5. 

In  all  cases,  only  satellites  which  are  known  to  contribute  to  the  FICA  file  are  allowed  to  be 
referenced  in  the  STARPREP  runstream.  If  a  particular  satellite  is  listed  in  the  STARPREP  runstream  and 
this  satellite  is  not  present  in  the  FICA  file  being  processed,  the  preprocessor  will  fail.  The  overall 
STARPREP  processing  time  for  a  typical  FICA  flle  is  only  about  2  to  5  minutes  (on  Unisys  equipment). 


STARPRCr  PLATTVIULE  ORT  £23  t**’ 


SHOC.P  •«LKCLASari£0** 

SMDG.P  *  *LNCLAS SI  PIED** 

9PRT  .S 

ESGT*HALYS.NCHSTAR/RUNbALL 

8PRT.S 

E$CT*MM.TS.MEVSTAR/RUMPALL 

iXQT  CPS  ‘rtCASTMPREP.PTSTARPREP/PCLOCK 

EXOT  CPS  »riCASTMPRCP.PTSTARPREP/PCLOC 

ASCriLES 

L 

GE0STAR*STATICKL06. 

ASGPILES 

4 

6 1 OSTAR  •$  TA  nON  LOG  . 

Assrtics 

7 

Ge0STAR«GASPtAA«223. 

ASGPILES 

7 

GEOSTAR  *64521  AA9223> 

asgpiles 

10 

6E0STAR*PTPLA92236. 

ASGPILES 

10 

6E0STAR*PTPLA9223P. 

ASGPILES 

11 

GC0STAR*038SPLA9223T  . 

ASGPILES 

11 

6E0STAR*0385RlA92231  . 

ASGriLCS 

u 

GC0STAR«a68SPLA9223T  . 

ASGPILES 

12 

6e0STAR*068SDLA9223T  . 

ASCriLES 

13 

6E0STAR*a98SPLA9«23T  . 

ASGPILES 

13 

6E0STAR*C98  5PLA9223t  . 

ASGPILES 

14 

GC0$TAR*118  5PLA9  2E3T  • 

ASGPILES 

14 

6e0$TAR*1l85PU9223T. 

asgpiles 

15 

6 E0STAR*1 28  SPLA9  223r  . 

asgpiles 

15 

GCOSTAR*1285RLAf2c3T  . 

asgpiles 

16 

GeOSTAP*t38SPLA9  223T  . 

ASGPILES 

16 

6L0$TAR*138SPLA9223T  . 

ASCFIlES 

31 

G  COSrAR«H  T8  5PLA  9  223P  * 

asgpiles 

21 

6PS*EP0 389218. 

asgpiles 

32 

GE0STAR*RC8SPLA92E3P> 

asgpiles 

22 

GPS*EP06892  18. 

asgpiles 

33 

6  EOSTAR'S  V8  SPLA  9 223P  . 

ASGPILES 

23 

6PS*EPC9892lB. 

asgpiles 

34 

GeasrAR«S78  5PLA9223P  • 

ASCFILES 

24 

6PS*CF11B9218. 

ASGf  IlES 

41 

GE0STAR*HNET. 

asgpiles 

25 

6PS*CF12892 18. 

asgpiles 

42 

G COSTAR •PRC  V. 

ASGPILES 

26 

6P$*EP1389218. 

asgpiles 

43 

GEOSTAR  •PSA  T. 

ASGPILES 

31 

6eOSTAR*nT8SPLA9223P  . 

ASGFIICS 

44 

GE0STAR*PSTA« 

ASGPILES 

32 

G  COST4R*R  C8  5PLA9  223P  . 

asgpiles 

51 

C  EOS  TAR  •P  TO  3223 L  CUT. 

ASGPILES 

33 

GEOSTAR *8X8  SPU9223P  . 

asgpiles 

52 

GC0STAR*PT06223L0UT. 

ASGPILES 

34 

GC0STAR*ST8  5PLA9223P  . 

ASGPILES 

S3 

C  COSTAR  •PT09223L0Ur. 

ASGPILES 

40 

GP5*PC89218. 

asgpiles 

54 

GC0STAR*PT1  1223LeUT. 

ASCPILU 

41 

CEOSTAR*nnCT. 

ASGF2LES 

53 

CE0SrAR*PT12232  L0UT. 

ASGPILES 

42 

6E0STAR*RRCV. 

ASGPILES 

54 

6C0STAR*PT1 3223L0UT. 

ASGPILES 

43 

6E0STAR*MSAT. 

PPROCSEG 

tOlT 

TTAGCOR  DAT  AC  OR 

ASGPILES 

44 

GCOSTAR^MSTA. 

CNPUTCCa 

ID 

2R  ER  TC  SC  GR  SA  TH 

ASGPILES 

5l 

GE0STAR*PT03223l0UT. 

APPLTCC* 

ID 

JR  ER  TH  SC  GR  SA 

ASGPILES 

52 

GEOSTAR*PT06223  l0UT. 

PLOTCCRS 

ID 

IR  ER  TH  SC  GR  SA 
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Figure  4.  Broadcast  and  Precise  STARPREP  Runstreams 
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Figure  S.  STARPREP  Runsoreams  for  Excluding  a  Satellite  (PRN  9  has  been  excluded) 

6  J  Broadcast  or  Precise  Ephemerides  and  Satellite  Clock  States 

The  STARPREP  preprocessor  allows  the  user  to  choose  from  two  sources  of  GPS  ephemerides 
and  satellite  clock  state  estimates.  These  two  sources  are  referred  to  as  broadcast  and  precise.  Thebioadcast 
ephemerides  and  clock  states  ate  predictions  which  are  derived  from  the  S-siaticm  network  operated  by  the 
US  Air  Force  while  the  precise  ephemerides  and  clock  states  are  weekly  nost-fit  estimates  computed  at 
DMAHTC  using  the  ON^S  inogram.  The  precise  ephemerides  and  clo^  states  are  estiinated  fiom  data 
collected  at  the  5  DMA  fixed  sites  and  the  5  Air  Force  sites.  This  combined  10-station  network  provides  a 
globally  distributed,  balanced  stadon  set  and  thus  considerably  improves  the  estiinated  ephemerides  and 
clock  states.  While  the  broadcast  and  the  precise  ephemerides/clocis  are  both  generated  using  the  WGS  84 
system  of  constants  and  models,  the  geodetic  user  must  remember  that  the  broadcast  messages  are 
[vedicdons  of  where  the  satellite  will  be  in  the  future  and  can  never  be  as  precise  or  accurate  as  the  post-fit 
estimates  which  are  generated  by  OMNIS.  The  disadvantage  of  using  the  precise  qihemeiides/clocks  is  the 
delay  of  a  few  weeks  before  the  precise  ephemerides/clocks  are  available.  This  is  rarely  a  problem  in 
geodedc  surveying. 

The  STARPREP  algorithm  is  designed  to  use  precise  satellite  clock  states  which  were  esumated 
along  with  the  precise  ephemerides.  The  only  time  the  broadcast  clock  states  are  used  is  when  the  broadcast 
ephemerides  are  selectecL  This  selection  occurs  automadcally  in  the  STARPREP  algcuithm. 

One  of  the  reasons  for  maintaining  the  option  to  process  broadcast  ephemerides/clocks  is  the 


possibility  of  using  STARPREP/GASP  in  the  field  on  an  appropriate  portable  computer.  Ineld  use  would 
require  the  use  of  broadcast  ephemerides/clocks.  While  the  positions  estimated  in  the  field  would  be  inferitH' 
to  those  estimated  from  the  precise  ephemerides,  they  would  be  adequate  for  validating  the  data  set  and 
would  provide  initial  coordinates  for  the  fixed  end  of  a  baseline.  In  the  usual  post-processing  environment, 
the  use  of  broadcast  orbits  can  facilitate  the  resolution  of  problems  with  a  data  set  and  provides  an 
opportuni^  to  evaluate  the  broadcast  data  if  the  station  coordinates  are  precisely  known  in  advance.  For 
routinct  precise  geodetic  point  positioning,  the  DMA  precise  ephemerides  and  satellite 
clock  states  must  be  used. 

6.4  The  Point  File 

The  primary  file  which  is  generated  by  STARPREP  is  referred  to  as  a  POINT  FILE’.  This  point 
file  contains  the  corrected  Li  data,  the  interpolated  ephemerides  and  satellite  clock  states  and  identifying 
information  which  indicate  the  corrections  that  were  applied  and  the  source  of  the  ephemerides  and  clock 
states.  A  sample  of  the  ASCII  form  of  a  point  file  (generated  by  the  utility  program  PTDUMP)  is  given  in 
Figure  6.  Note  that  the  point  file  form  used  by  STARPREP  is  binary.  The  structure  for  the  contents  of  a 
point  file  is  also  given  in  Figure  7.  The  ASCII  form  can  be  visually  inspected  if  problems  occur. 

A  POINT  file  naming  convention  has  been  established  within  DMA.  The  structure  of  the  name  is 
as  follows: 

Example:  GEOSTAR*PT2719148PA. 

Where  'GEOSTAR*  is  the  Unisys  qualifier  and  would  not  be  necessary  on  non-Unisys  equipment. 

The  remainder  of  the  file  name  is  die  critical  part' 

Character  1-2  'PT  indicates  a  POINT  file 

3-S  '271'  last  3  digits  of  station  identifier 

6  *9'  last  digit  of  year  data  was  ctdiected 

7-9  '148'  day  number  of  year  that  data  was  collected 

10-12  PA'  T  indicates  Precise  q>hemerisA:locks. 

'B'  indicates  broadcast  ejAemeris/clocks 
'A'  (optional)  or  some  other  designator 
indicates  which  subset  of  the  original  data  set 
was  used 

For  DEC  (VAX  or  MicroVAX,  the  file  name  extension  'PNT  should  be  used. 


.cctceaoo  C  «.«S74t«V<T*!3t71*« 

t  7  71  »J 

13  S«  «S  !3 

ll•lll«1Csc»M*7*  lo.acooiie 

•  *1«  ,',«71t.3i3«1732''lf3 

•:77i.i2C7i«iM(94f*  •a2a«.ti;e3:i««!*it 

-.1C<(S*5T27271I«I7»-Ct2  - .ZS472I3C11  CIlllCU-SOl 

.l*1J«443!in«f  .S404*73(l31Z!S!:iC-«Q2 

•  413  11*74«.ei!«47342324 

-!74  <.113S44«7  ii«M4  -2o2TZ.44Cla«CTtC12t 

..1Ct3411f|3l«4iai1S>C02  -.3SZSa71<«3(lSI<7!«-0C2 

•.3S  iini2!33«ii3:4«>ca2  Jivzoitiivitiviice-qoa 

4  41«  1la73«.Cl311«4414SS 

.ZT3:.7Z47':44Sa«S132  •■103SS.«S(4«44239«I1 

-.iel<4SS<l3eS44Z<01*{02  - J4«S44S4t4|43lfS44-Ca2 
.!SI!S3t3f1U7S4iat>«0Z  JSr7SSa043 14104437-00: 

4  41*  11tlZ«.C13 7**01*370 

-174  7.1*«Z;3«*330IS2  -a3*0.«l4f1 1  lltlllz 

1010471047 T4*77SI*'aoZ  -.*la23«71*7!1<l!*11>00Z 
-.31  r|C]21t41«*73  •34-002  d3*31f*e2*  i*«oi*«o-aoz 

77  *7*  ,taS3*.lC<74<14C1]* 

-12074.4112147104441  -203*7.13411 !>4031*7 

-  .11  •3*703*707174473-004  -  .33**3*112a*133«37*-a0Z 
^44-4*S*l<144477724-ce2  .441e7*f1|Ti4l*4;i«-002 

•  41*  111***. 013*14303*17 

.l7!l.t**00*!>03**!13  -  204*7.4*2*02*134333 

-.17 144771  iZ7z*73ti*-eei  -  a  loaiiiis  itct  3*12-002 

.31  ••*4113723123  :lC-oel  ^3**2:3720437101*1-002 


210.441* •*•**0**33* 
13 
ft 


7**00.40  7*  2*  *1471 7» 
1442C.7S011 34*44147 
.4  071* 073 44  40 tt  30  3  3 -O  0  3 
.3000000000 COCOOBOO 
1*113.  332:71174244* 


1*3*3.3910*7021  71*1 
.403702*22307***113-003 
.00000000 00 00 000000 
1*1 TO .3*37 7147*20*1 
**4«*.*3010**31C*0* 
.4C33***043 1**4 37 1 *-003 

jiaocoaoaooooooeooo 

1*133. 12 40*3*1*2311 
1*5**.03* 31 5*2172** 
.311114171S  23**02  *0-403 
.400  cocao 00 00900000 
21031. 11a>  471112334 
12011.710  3*077*3401 
.20132473 2* 11412412-0 93 

.ocacoBoowceoeceoo 

1«149.n474134*4«l3 

1*117.43020*79*7371 

.3**c2lSSl3»*42»71-0eS 

•ocococnaococccooo 


.:i*««*«***«*<***t*-oai  *  1 

0  C  2  7  II 

00*27  114*  7  7llIt(4TMtC**I 

-2*3  .c313333*S«27l*0  .OOOOOOOOaoCaCDOOOO 

.490  rcccao'joacaooea  -.112**34**12*4 11*24.001 

.272l<7333t**14*30e-Ca3  -141.ai4«40***7t770* 
.2it<4:3i-co7  ,i23«oai*-eii 
-23I.C*S2*«77704**37  30.0000000000000000 

.000  ccecoo  eoocooooo  -.11303«2**72140«*  07-001 

.*727741730***13*21  -003  -24X.0iril7721*M17* 
.131103*7-007  .ii*o2*a*-aii 
-  273  .0 43*1 121  *14*1 1 1  *0.0000000400000000 

.nococcooeoe  090009  -.7i3ets3o*4iii4*407-on 

.I710IC**1**«**1707-OOS  -143.0172744*1*0331* 
.173771*1-007  .32*330**-011 

-Z17.»l2ei«377Z*l270  *0.0000000900000000 

.000000000000090001  -.113134344***33o3*l-«91 

.«731t7a3433  *  0»0n  1  -003  -243.01743172*4 711«4 
.1Z3*722*-9o7  .321|e3r9-01i 

-24. *0**017*133*3*31  .100000000100000000 

.oooccccooeaacoooQo  -.247*73*1*171*  3*23*-«oi 

.«**1'C*20*42I1**21  -003  .24.3*e72t7**S17**7* 
.13*02132-007  .  329(c**4-e11 

-3O2  .*I3*37037*«2»27  120.000000000000000 

.000900000000000000  ..11311377131137*111-001 

.173  2*32***39444443  -003  -243.ft171***0**1*c2* 
.244l041*-ca7  .3217*343-112 


Figure  6.  Sample  of  a  Portion  of  a  Point  File 
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PCIM  flLE  RECORl  DESCRIPTIOK 


RECORt  »1  (riLf  header.  1  PER  fILcJ 


UOfiO  9 

-  VARI/BLe  - 

Trpc  - 

CESCrIpTION 
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STI  t 
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STATION  MIAEeR 
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SIT  ID 
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STATION  CCCRDINATf,  LATITUDE  tOEG) 
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STATION  CCtRDINATE,  LOK&ITUDE  fOEC) 
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LIST  or  SATELLITES  IN  DATA  FILE 
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SVR  ECN 
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ACDRESS  FCF  EACH  SATELLITE'S  DATA  fAOT  CSEC) 

23 

CAT  tP 

C«6 

PROCESSING  DATE 

2A 

ITR 
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YEAR  OF  START  OF  DATA 

25 

HSTTPE 

C*2 

NEASUREHEAT  TYPES  ON  F IL c 

26 

EPHTTP 

C«2 

EFNEHERIS  TYPE 

27 

NCORRA 
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NtPEER  OF  CORRECTIONS  APPLIED 

2S 

CCRFA 

C*2 

LIST  OF  CCRRECTIONS  APPLIED 
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NIMEER  OF  CORRECTIONS  COMPUTED 

30 
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C,2 

LIST  OF  CCRRECTIONS  COPPUTED 

3t 

RRN11 

S 

REDLESTED  PEASURERENT  INTERVAL,  SeC 

52-5; 

COHPNT 

C*3c 

CCHPENTS  CN  FILE  HEADER 

RECCRt  *2  (DATA  RECORDS,  1  FOR  EvErT  CATA  bcImT) 


UORD  f 

-  VARI  ABLE  “ 

TYPE 

-  CESCRIPTION 
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PRN 
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SATELLITE  FRN  NUMBER 
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CCRKCCTED  COPPLER  HEASUREKENT, 

i 

TcUP 

0 

CtPUlATlVE  TIPE 

7 

XSAT(T) 
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XSaT(3) 

f» 

lATERPOLATFC  SATELLITE  POSITION.  I-COP.POAEHT 

tc 

OOPCFF 

c 

OCPFLER  OFFSET  CORRECTION,  KH  (NOT  USED) 

IT 

EROT 

e 

Earth  rotation  correction,  kp 
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D 

GENERAL  RtlATIVlTY  CORRECTION,  KP 
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21 
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DCPFLER  PtASUREPENT  VARIANCE,  KP.«*2 

Figure  7.  Point  File  format 
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6.5  Quality  Control  -  Analysis  of  STARPREP  Output 

Selected  output  from  a  normal  STARPREP  run  is  shown  in  Appendix  B.  The  collection  of  input 
flies  which  were  used  to  generate  the  output  ^own  in  Appendix  B  is  available  as  a  test  data  set.  The  input 
flies  consist  of  the  FICA  file,  the  DMA  precise  ephemeris  flies,  and  the  precise  satellite  clock  file.  The 
PLOTCORR  option  should  always  be  invoked  in  the  STARPREP  runstream.  This  option  produces  plots 
of  all  applied  data  corrections,  by  satellite,  as  part  of  the  STARPREP  ouqrut.  These  plots  are  very  helpful 
in  understanding  the  corrections.  They  are  also  very  helpful  in  resolving  problems  with  a  particularly 
difficult  data  set  Note  that  the  plots  do  not  contain  values  for  every  data  point  Only  every  nth  point  is 
used  to  generate  the  plots.  Each  STARPREP  output  plot  should  be  examined  for  the  following 
characteristics: 


to  check:  (Note  that  all  corrections  plotted  in  the  output  are  in  units  of  kilometers) 

>  TJncorrected  range  vs  dme'  and  TJncorrected  Doppler  vs  time':  The  plots  should  be  linear  or  2nd  order 
curves  with  possible  gaps  where  no  data  was  collected  from  a  particular  satellite.  Discontinuities  in  the 
Doppler  plot  should  correspond  to  cycle  slips  or  scenario  changes. 

>  ‘Satellite  #  range  iono  vs  time':  Corrections  plotted  on  this  graph  should  range  between  0  and  -IS 
meters.  All  corrections  should  be  negative.  The  random  nature  of  these  corrections  reflects  the  normal 
pseudotange  noise. 

>  ‘Satellite  #  Doppler  iono  vs  time':  These  ionospheric  corrections  are  used  to  identify  cycle  slips  in 
either  Li  or  L2.  These  cycle  slips  appear  as  steps  in  the  correction  values.  While  these  cycle  slips 
normally  do  no  harm  in  the  GASP  algorithm,  they  ate  of  interest  for  further  detailed  analysis.  These 
corrections  should  be  in  the  range  of  0  to  ±1S  meters. 

>  'Satellite  #  sat  clock  vs  time':  This  graph  should  appear  linear  or  vary  smoothly  over  the  span  of 
data.  These  corrections  are  interpolated  from  the  precise  satellite  clock  biases  and  serve  to  correct  the 
individual  satellites  to  GPS  time.  When  txoadcast  ephemerides  and  clocks  are  used,  these  corrections 
ate  computed  flom  the  polynomial  coefficients  found  in  the  FICA  Block  9s. 

>  'Satellite  #  relativity  vs  time':  This  plot  should  vary  smoothly  over  the  data  span.  The  plot  should  be 
linear,  or  a  2nd  or  3rd  order  curve.The  correction  values  should  range  between  0  and  ±15  meters. 

>  'Satellite  #  earth  rotation  vs  time':  This  plot  should  be  linear,  or  a  2nd  or  3rd  order  curve.  The 
correction  values  should  be  in  the  range  0  to  ±40  meters. 

>  'Satellite  #  sat  antenna  v$  time':  This  small  correction  should  be  linear,  or  a  2nd  or  3td  order  curve. 
The  correction  values  should  be  less  than  1  meter. 

>  'Satellite  #  Hopfleld  troposphere  vs  time':  The  plot  should  be  linear  or  a  2nd  or  3rd  order  curve.  This 
correction  is  dominated  by  the  elevation  angle  and  as  such,  should  vary  smoothly  over  the  span.  A 
spike  or  discontinuity  may  indicate  that  an  incorrect  ephemeris  file  was  assign^  in  the  runstream. 
Correction  values  should  be  in  the  range  0  to  30  meters. 

Othg  items  lo  chwit: 

>  Check  that  the  Ephemeris  Source'  is  Reference'  if  the  DMA  precise  ephemerides  and  clock  states  were 
selected  in  the  runstream.  In  the  rare  cases  when  STARPREP  encounters  a  problem  locating  the  proper 
precise  q^hemeris  and  satellite  clock  flies,  it  will  default  to  the  Broadcast  ephemeris  and  clocks.  This 
will  be  noted  in  the  Ephemeris  Source'  message. 

>  Check  that  the  'Corrections :  Computed  Applied'  are  as  expected. 

>  Check  and  note  'Output  File  Name'  as  the  Point  File  name  which  will  be  used  to  input  to  GASP. 
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7.  POINT  POSITION  ESTIMATION  (GASP) 

7.1  Overview  and  Sample  Runstream 

Two  files  which  STARPREP  produces  are  required  inputs  to  the  point  positioning  program 
(GASP)  (The  Qeodetic  Absolute  S.equential  Eositioning  Program).  These  files  are  the  point  file  and  the 
station  file.  The  only  information  needed  from  the  station  file  is  the  antenna  height  (monument  to  base 
of  antenna).  STARPREP  obtains  this  antenna  height  from  the  FICA  file  Block  101  or  Block  124.  Recall 
that  the  antenna  height  should  be  verified  against  the  field  sheets  when  the  data  set  is  reviewed.  GASP 
automatically  applies  the  constant,  known  offset  from  the  base  of  the  TI  antenna  to  the  mean  electrical 
center  (0.22733  meters). 

The  standard  GASP  runstream  is  much  simpler  than  STARPREFs.  Only  two  files  are  assigned 
and  unless  some  experiments  are  being  conducted,  the  standard  options  (A)  are  all  that  need  to  be  specified. 
An  example  of  the  standard  GASP  runstream  and  an  example  of  a  GASP  runstream  which  allows  selection 
of  non-standard  options  are  shown  in  Figure  8.  A  complete  sample  GASP  ouqwt  is  given  in  Appendix  C. 

Typical  overall  processing  time  for  GASP  is  only  1  to  2  minutes  (Unisys  equipment). 
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Figure  8.  Sample  GASP  Runstreams 
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7.2  Algorithm  Description 


The  corrected  carrier  beat  phase  data  is  known  to  contain  biases,  integer  cycle  ambiguities  and  other 
undesirable  characteristics  which  are  dependent  on  factors  such  as  initial  acquisition  epoch,  cycle  slips  and 
receiver  frequency  standard  fluctuations.  Moreover,  common  errors  may  be  iriadvertently  introduced  through 
the  application  of  less  than  perfect  data  corrections  such  as  those  associated  with  the  atmosphere,  and  less 
than  perfect  estimates  of  the  satellite  ephemerides  and  clock  states.  In  an  attempt  to  remove  these 
undesirable  contributions,  the  GASP  algorithm  is  designed  around  a  differencing  scheme  which  cancels 
common  errors  in  the  corrected  data  and  adapied  satellite  ephemerides  and  clock  states. 

The  first  step  in  forming  an  observable  is  to  difference  two  consecutive  carrier  beat  phase 
observations  (convert^  to  kilometers)  of  the  same  satellite.  Typically,  the  interval  between  consecutive 
observations  is  30  seconds.  This  between-epoch  difference  (a  bia^  delta-range)  is  then  differenced  with  a 
corresponding  between-epoch  difference  from  another  satellite.  One  'GASP  observable'  is  formed  from  data 
collected  at  one  station  tom  two  satellites  at  two  consecutive  epochs.  Since  the  TI4100  receiver  tracks  up 
to  four  satellites  simultaneously,  a  given  pair  of  data  epochs  can  yield  up  to  three  'GASP  observables'.  A 
graphic  representation  of  this  differencing  scheme  is  given  in  Figure  9.  For  each  pair  of  data  epochs 
utilized,  one  satellite  is  used  as  the  reference  tom  which  the  others  are  differenced. 


Figure  9.  Formation  of  GASP  Observables  (A  &  B  represent  two  epochs  separated  by  a  30  second  interval) 

After  forming  an  array  of  GASP  observables,  a  least  squares  estimation  technique  is  applied  such 
that  the  vector  of  estimated  parameters  contains  an  'alteration'  to  the  a  priori  Earth-centered,  Earth-fixed 
Cartesian  station  antenna  position  components.  No  other  parameters  are  estimated.  After  three  iterations  of 
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the  non-linear  model,  the  estimated  parameters,  their  scaled  variance-covariance  matrix  and  the  RMS  of 
residuals  are  passed  to  a  sequential  estimation  algorithm.  This  sequential  technique  (a  Kalman  filter) 
utilizes  the  RMS  of  residuals  from  the  least  squares  fit  as  the  variance  of  a  GASP  observable.  In  the 
sequential  estimation  technique,  a  parameter  update  is  performed  after  each  observable  is  processed.  This 
allows  plots  to  be  generated  which  show  the  position  component  estimates  as  a  function  of  the  data  set 
span.  The  level  of  convergence  in  the  plots  is  one  indication  of  the  precision  of  the  position  estimate. 
Analysis  of  residual  plots  and  of  a  posteriori  variance-covariance  matrices  are  also  helpful  in  evaluating  an 
individual  result. 


7.3  Pseudorange  Minus  Carrier  Beat  Phase  Biases 

Before  the  GASP  observables  are  formed,  each  satellite  is  analyzed  for  consistency  between  its 
pseudorange  data  and  its  carrier  beat  phase  data.  A  delta-range  is  comput^  for  each  pair  of  data  epochs  by 
subtracting  two  corrected  pseudorange  observations  and  also  by  subtracting  two  corrected,  convert^  carrier 
beat  phase  observations.  Theoretically,  these  two  delta-ranges  should  be  equal.  Since  however,  there  are 
different  noise  levels  associated  with  these  two  data  types,  they  will  not,  in  general,  be  equal.  A  mean 
difference  between  these  two  types  of  delta  ranges  is  computed  for  each  satellite  in  a  point  file.  These  mean 
differences  are  listed  in  the  beginning  of  the  GASP  output  These  mean  differences  are  referred  to  as  delta 
pseudorange  minus  delta  carrier  beat  phase  biases.  These  mean  values  are  removed  from  each  observable 
tested  in  the  procedure  which  edits  data  based  on  a  5  meter  tolerance.  (They  ate  removed  only  fm-  the 
tolerance  test  not  in  the  construction  of  the  observable.)  Occasionally,  the  list  of  these  biases  will  reflect  a 
problem  with  a  particular  satellite  in  the  data  set  Quality  control  criteria  regarding  these  biases  is  presented 
in  Section  7.8. 

7.4  Reference  Satellite  Sequencing 

Whenever  the  GASP  standard  options  are  used,  the  reference  satellite  is  selected  sequentially,  such 
that  for  every  new  pair  of  epochs  proces^,  the  next  higher  available  PRN  number  is  used  as  the  reference 
in  the  differencing  scheme.  The  choice  cycles  back  to  the  lowest  available  PRN  number  after  the  list  of 
tracked  satellites  has  been  exhausted. 

Before  adopting  this  sequential  reference  satellite  sdection  process,  the  authors  tested  the  concept 
that  the  satellite  with  the  most  stable  clock  should  be  used  as  the  reference  throughout  the  data  span. 
Comparison  of  statistics  from  repeated  estimates  however,  indicated  that  the  sequential  approach  offered 
many  more  benefits.  These  benefits  are  rationalized  in  terms  of  reduced  correlation  among  the  observables. 
The  implementation  of  this  sequential  referoKe  satellite  selection  process  improved  day  to  day  repeatability, 
the  RMS  of  residuals  (for  most  fits),  and  the  variance-covariance  and  correlation  matrices  of  the  estimated 
parameters  by  up  to  30%.  In  particular,  the  G  ASP-estimated  longitudinal  component  reaped  almost  all  of 
the  beneflts.  There  remains  an  option  in  GASP  which  allows  the  user  to  override  the  automatic  sequencing 
of  reference  satellites.  This  option  is  only  meant  for  experimental  purposes. 

7.5  Rejected  Data 

Two  kinds  of  data  editing  are  performed  by  GASP.  The  first  kind,  discussed  briefly  in  Section  7.3, 
relies  on  a  comparison  of  delta  pseudoranges  with  delta  'converted  carrier  beat  phase'.  If  these  two  kinds  of 
delta  ranges  differ  by  a  user-selected  tolerance  (S  meters),  the  GASP  observable  (composed  of  two  delta 
converted  carrier  beat  phase  observations)  is  rejected  from  the  data  seL  This  kind  of  editing  is  referred  to  as 
Dseudorange  editing.  A  cecond  editing  method  uses  an  RMS  screen  of  3  times  the  previous  iteration's  RMS 
of  residuals  (the  RMS  is  initialized  before  the  first  iteration).  Any  observables  which  have  a  residual  greater 
than  3  times  this  RMS  are  rejected  from  the  data  set. 

Experience  with  GASP  has  shown  that  the  pseudorange  editing  process  generally  performs  most  of 
the  editing.  The  RMS  screening  process  usually  only  edits  significant  amounts  of  data  when  the  a  priori 
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coordinates  are  very  poor  or  if  a  particular  satellite  is  behaving  erratically.  Quality  control  criteria  regarding 
rejected  data  are  given  in  Section  7.8. 


7.6  Batch  Least  Squares  versus  Sequential  Estimation 

As  described  in  Section  7.2,  GASP  utilizes  two  different  esdrnation  methods.  The  first  is  a  'batch' 
least-squares  technique  and  the  second  is  a  sequential  technique,  usually  referred  to  as  a  Kalman  filter.  Since 
each  of  these  method  offers  specific  advantages,  they  are  used  together  to  maximize  the  extraction  of  useful 
information  from  the  adjustment  process.  The  standard  options  in  GASP  employ  three  iterations  of  the 
batch  least  squares  method,  then  one  iteration  through  the  sequential  method.  Results  from  the  batch 
method  provide  all  the  inputs  to  the  sequential  method.  The  estimated  station  coordinates  from  the  batch 
method  serve  as  a  priori  coordinates  for  the  sequential  method.  The  uncertainty  on  the  coordinates  entering 
the  sequential  stage  are  de-weighted  to  avoid  constraining  them  to  the  batch  results.  The  standard  deviations 
are  multiplied  by  100  before  they  are  used  in  the  sequential  processor.  Most  importantly,  the  batch  method 
provides  a  measure  of  data  noise  to  the  sequential  processor.  The  final  RMS  of  residuals  from  the  batch 
processor  is  used  as  the  standard  deviation  of  one  GASP  observable  processed  through  the  sequential  filter. 
Any  data  that  was  rejected  in  the  batch  processor  remains  rejected  in  the  sequendal  processor. 

Since  the  sequential  processor  updates  the  estimate  statitm  coordinates  as  each  new  observable  is 
processed,  a  series  of  estimated  station  coordinates  can  be  plotted  as  a  fiincdon  of  the  data  span  used.  These 
plots  are  provided  as  part  of  the  standard  options.  The  level  of  convergence  in  these  plots  (one  each  for  the 
X,Y  and  Z  components)  serves  as  an  indio^n  of  the  precision  of  the  final  estimated  coordinates. 

7.7  Variance-Covariance  and  Correlation  Matrices 

The  batch  least  squares  and  sequential  estimation  techniques  both  provide  for  the  estimation  of 
uncertainties  on  the  estimated  parameters.  In  our  qiplication,  the  estimated  parameters  are  the  three 
Earth-centered,  Earth-fixed  (ECEF)  Cartesian  station  coordinates.  An  estimate  of  uncertainty  on  these 
coordinates  can  be  obtained  firom  the  variance-covariance  matrix  for  the  estimated  parameters.  A  measure  of 
linear  independence  among  the  estimated  parameters  is  obtained  by  developing  the  correlation  maoix  of 
estimated  parameters.  The  variance-covariance  matrix  and  the  correi^on  matrix  are  symmetric  matrices. 

A  posteriori  standard  deviations  for  the  estimated  ECEF  station  coordinates  are  obtained  by  taking 
the  square  root  of  each  diagonal  element  of  the  variance-covariance  matrix.  These  diagonal  elements 
represent  the  estimated  variance  of  each  ECEF  component  Experience  has  shown  that  these  standard 
deviations  are  realistic  estimates  of  the  precision  of  the  estimated  point  position  components.  The  user  will 
notice  very  good  agreement  between  the  standard  deviations  estimated  in  the  batch  least  squares  processor 
and  the  sequential  processor. 

The  off  diagonal  elements  of  the  variance-covariance  matrix  are  the  covariances  and  are  used  to 
compute  the  linear  correlation  coefficients  which  make  up  the  correlation  matrix. 

7.8  Quality  Control  •  Analysis  of  GASP  Output 

The  output  from  each  GASP  run  is  to  be  analyzed  to  ensure  that  the  estimated  position  and 
relevant  statistics  are  reliable.  The  following  list  describes  the  items  which  must  be  checked. 

>  The  selected  options  that  will  be  used  are':  Check  all  values.  For  normal,  standard  processing,  the 
reference  satellite  should  be  00  to  indicate  reference  satellite  sequencing. 

>  'Echo  of  station  file'  (unit  34)':  Check  station  file  name  and  a  priori  station  coordinates. 
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>  'Distance  From  Ground  Mark  to  Base  of  TI4100  Antenna':  Check  against  height  from  field  notes. 

>  'Summary  of  Reference  Satellite  Selections':  All  PRN's  in  the  associated  STARPREP  run  should  be 
shown.  When  reference  satellite  sequencing  is  used  (as  it  should  be  in  all  standard  runs)  the  "No.  of 
Times  Used'  should  be  approximately  commensurate  with  the  amount  of  data  available  from  each  PRN. 

>  'Mean  of  Delta  Pseudorange  -  Delta  Carrier  Beat  Phase  for  Each  Satellite  in  this  Data  Set':  In  general,  a 
good  data  set  will  yield  values  for  these  biases  which  are  in  the  range  of  a  few  centimeters  to  a  few  tens 
of  centimeters.  Values  higher  than  about  50  centimeters  indicate  a  problem  with  that  particular 
satellite  or  a  small  amount  of  data  from  that  particular  satellite.  If  any  of  these  biases  is  orders  of 
magnitude  different  from  those  of  the  other  satellites,  go  back  to  STARPREP  and  exclude  the  outlier 
satellite  from  the  data  set. 

>  Rejected  Data':  Whenever  an  observable  is  rejected  by  GASP  the  following  information  is  provided  to 
the  user. 

Reason  for  rejection:  PSR  EDn'=>  Pseudorange  edit  tolerance  failed;  RMS  EDrr=>  Residual  outside  of  RMS  screen 
Observable  location  in  the  dau  array  (DPAR)  (two  locations,  this  array  is  arranged  in  matched  pairs) 

Satellite  PRN  numbers  (two  different  PRN  numbers  representing  the  between  satellite  difference) 

Time  tags  (these  should  always  differ  by  the  data  collection  rate,  usually  30  seconds) 

Difference  of  delta  ranges  for  each  satellite  (for  a  PSR  edit) 

Residual  and  RMS  screen  limit  (for  an  RMS  edit) 


In  very  rare  cases,  large  systematic  values  for  a  PRN#  in  PSR  EDIT>  arc  found.  This  condition 
indicates  a  problem  with  the  data  (anomalous  pseudorange  or  carrier  beat  phase  data  values  have  entered  the 
Point  file).  Since  each  PRN  number  is  sequentially  paired  with  other  PRN  numbers  to  form  the 
observables,  a  large  amount  of  data  from  'good'  sateUiies  can  be  thrown  away.  If  the  anomalous  data  values 
cannot  be  located  and  removed,  the  satellite  must  be  excluded  in  the  STARPREP  run.  This  problem  is 
sometimes  eliminated  by  altering  the  data  span  by  30  seconds. 

7.8.1  After  Three  Batch  Iterations  Check  the  Following: 

>  'A  Posteriori  Variance  of  Unit  Weight':  Should  be  in  the  range  0.8  to  1.2.  This  insures  that  a  valid 
measure  of  data  uncertainty  has  been  passed  to  the  sequential  estimation  module. 

>  'Correlation  Matrix  of  Cartesian  Estimated  position':  If  any  off-diagonal  correlation  coefficient  is  near  1 
in  absolute  value  (>0.8),  a  problem  should  be  suspect^.  The  geometry  of  the  data  set  should  be 
analyzed  along  with  the  data  rejected  by  GASP. 

>  'Estimated  Standard  Deviations  on  X,Y,Z  (m)':  These  standard  deviations  should  be  less  than  1.5 
meters.  The  Z  component  is  usually  smaller  than  X  or  Y. 

>  'Current  RMS  of  (Mean(O-C)-(O-Q)':  Should  be  less  than  5  cm.  'OC  is  Observed  minus  Computed 
where  'O'  is  the  GASP  observable  constructed  from  the  corrected  LI  carrier  beat  phase  data  and  C  is 
the  predicted  (modeled  .computed)  value,  obtained  from  knowledge  of  the  satellite  state  vectors  and  the  a 
priori  station  position  vector. 

>  Number  of  GASP  observables  Used':  Should  be  greater  than  220. 

>  'Percentage  of  Data  Rejected':  After  each  iteration,  the  percentage  of  rejected  data  is  given.  Typically, 
this  percentage  is  between  about  5  and  15%.  Before  each  new  iteration  is  started,  a  summary  of  the 
remaining  data  epoch  pairs,  sorted  by  satellite,  is  listed.  By  analyzing  the  rejected  data  and  the 
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remaining  qxxJi  pairs  list,  it  may  be  found  that  a  particular  satellite  is  being  nearly  completely  rejected 
from  the  dam  sec  If  more  than  23%  of  the  data  was  rejected,  go  back  to  STARPR^  and  regenerate  the 
point  file  with  the  troubled  satellite  excluded. 

>  Check  the  Plot:  'OMC(km)  vs  Time  (secs)  for  #  3  Iteration'  These  'residuals'  must  be  of  a  random 
nature  on  the  3rd  iteration.  Note  that  in  the  first  iteration,  poor  a  priori  station  coordinates  will 
produce  systematic  traces  on  these  plots,  this  condition  presents  no  problem  as  long  as  a  normal 
amount  of  data  was  rejected.  The  'O-C  plot  on  the  third  iteration  must  always  appear  random.  Bad 
data  or  incorrect  ephemoides  will  also  produce  systematic  traces. 

7.8.2  After  the  Sequential  Filter  Processing,  Check  the  Following: 

>  All  items  listed  in  section  7.8.1  (there  will  only  be  one  'OMC  vs  Time'  plot  and  no  a  posteriori 
variatKe  of  unit  weight) 

>  "Mean  of  Delta  pseudorange  -  Delta  carrier  beat  phase  for  each  satellite  in  this  data  set':  as  described  in 
section  7.8. 

>  Plots  of  'X,Y  and  Z  Coordinate  (km)  vs  Hme  (secs)':  These  important  plots  show  how  the  Kalman 
filter  estimate  converges  over  the  data  span.  If  any  of  these  components  appear  to  'wander  around' 
beyond  the  level  of  1  meter  after  a  few  hours  of  data  has  been  processed,  the  user  should  suspect  a 
problem  and  begin  examining  the  data  set  The  user  might  find,  for  example,  that  the  receiver  tracked 
only  2  satellites  ftu'  some  period  of  time.  This  period  may  correspond  to  the  period  where  the  station 
coordinate  estimates  are  'wandering'.  If  such  a  situation  occurs,  the  data  set  should  be  edited  (to  keep 
only  periods  where  4  or  more  satellites  were  tracked)  or  discarded. 

>  Compare  the  filter's  position  estimate  to  the  estimate  given  at  the  end  of  the  batch  processor.  The 
numeric  comparison  is  already  done  where  the  'Initial  •  Final  Station  Differences'  are  given.  In  the  case 
when  the  sequential  processor  is  used,  the  'initial'  station  position  estimate  is  provided  by  the  batch 
least  squares  processor.  These  two  position  estimates  should  not  differ  in  any  component  by  more  than 
30  cm.  Theoretically,  if  no  further  ^ta  is  rejected  in  the  sequential  processor,  and  if  the  'Q'ing  is  set  to 
zero  in  the  Alter,  th^  two  position  estimates  should  be  eqi^. 
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8.  PRECISION  AND  ACCURACY 

8.1  Deflnitions 

The  terms  'precision'  and  'accuracy'  have  distinct  meanings  and  should  not  be  interchanged  in  any 
discussion  of  geodetic  positioning  results.  The  precision  of  an  estimated  parameter  is  a  measure  of  the 
degree  to  which  the  results  are  repeatable.  In  general,  it  is  acceptable  to  interchange  the  terms 
'repeatability',  'intenud  consistency',  and  precision.  The  standard  deviations  of  the  estimated  station 
coordinates  which  GASP  generates  are  measures  of  precision. 

In  order  to  evaluate  accuracy,  a  well-kix>wn,  independent  standard  of  known  higher  accuracy  must 
be  available.  In  other  words,  some 'truth' values  must  be  available  for  comparison.  In  geodetic  positioning 
work,  these  'truth  coordinates'  might  be  obtairted  from  other  advatKed  satellite  geodesy  methods  such  as 
Satellite  Laser  Ranging  (SLR).  If  no  such  independent  standards  of  accq)ted  higher  accuracy  are  available, 
there  is  no  direct  way  to  evaluate  the  accuracy  of  positioning  results.  Even  when  highly  accurate  standards 
are  available,  reference  frame  differences  may  be  present  which  will  complicate  the  evaluation. 

Figure  10  gra^cally  portrays  the  distinction  between  accuracy  and  precision  ftv  a  2-dimensional 
example.  In  the  diagram  shown  on  the  right  side  of  Figure  10,  the  offset  between  the  cluster  of  estimates 
and  the  'truth'  position  is  referred  to  as  a  'bias'. 
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Figure  10.  Graphic  Example  of  Accuracy  and  Precision 

To  quantify  the  inecision  and  accuracy  of  GASP  position  estimates,  an  extensive  study  was 
conducted  by  Malys  and  Jensen  [1989].  Using  over  200  data  sets  and  the  standard  qrerating  proc^ures 
outlined  in  this  document,  the  overall  level  of  (xecision,  measured  by  the  standard  deviation  of  repeated 
positioning  results,  was  found  to  be  73  centimeters  in  each  component  An  accuracy  assessment  was 
performed  by  comparing  a  small  number  of  GASP  results  to  positioning  results  obtained  through 
VLBI/SLR  methods.  The  mean  difference  between  a  GASP  component  and  a  VLBI/SLR  component  in 
absolute  value,  was  76  centimeters.  These  levels  of  precision  and  accuracy  are  achievable  on  a  routine  basis 
by  following  the  guidelines  presented  here. 

8.2  Reporting  Positioning  Results 


Since  there  are  a  number  of  methods  available  for  geodetic  positioning  with  GPS  data,  results 
must  always  be  carefully  documented  and  cataloged  in  a  way  which  cle^y  shows  which  kind  of  positioning 
was  performed  (absolute,  static  relative,  or  kinematic  relative)  and  which  software  packages  were  used. 
Since  many  GPS  positioning  algorithms  are  still  under  development  and  undergoing  enhancements,  the 
software  version  numbers  should  also  be  recorded.  For  point  positions  estimated  by  GASP,  the  following 
minimum  specifications  should  be  entered  into  the  user's  data  base  or  whenever  the  results  are  provided  to 
another  organization: 
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Station  Name: 

Method: 

Software: 

Epbemerides  and  Satellite  Clock  States: 
Date(s)  of  Occupation: 

GPS  Receiver: 

Date  of  Point  Position  Estimation: 

Data  Collection  Span: 

PRN  Numbers  Tracked: 

Fiual  Number  of  GASP  Observables: 
Percentage  of  Data  Rejected: 

RMS  of  Residuals: 

WGS  M  Esrtmafgd  Station  CnnrHinates 

Monument  (X,Y^)  (meters): 
Standard  Deviations  (X,Y^)  (meters): 


DMA123 

GPS  WGS  84  Absolute  Point  Posititming 
STARPREP  (veision  1.0),  GASP(vcrsion  1.0) 
Precise  WGS  84  EF##89246.  PC89246 
Day  247, 1989 

TI4100  using  GESAR  (version  1.9) 

OS,  December,  1989 
4.5  Hours 
3,6,9.11,12,13 
350 
112% 

4.1  cm 


327259.205  -6340165.804  612494.684 

0.732  0.643  0.451 


Geodetic  coordinates  (monument):  ^  5®  32’ 51284"  X»  272°  5T  17.293"  h=  144.823  meters 

Antenna  Height  (monument  to  electrical  center  of  antenna):  1.523  meters 

The  standard  deviations  listed  above  can  be  used  in  number  of  ways  to  express  statistical 
probabilities  regarding  the  positioning  result  Under  the  assumption  that  the  txion  in  a  positioning  result 
are  normally  distributed,  one  standard  deviation  Otl*’)  rqpresents  68.27%  of  the  area  under  the  normal 
distribution  curve.  If  one  wishes  to  express  a  90%  confidence  level  for  a  positioning  result  the  standard 
deviations  listed  above  should  be  multiplied  by  1.6449.  Ninety  percent  of  the  area  under  the  normal 
distribution  curve  is  represented  by  1.6449  a.  For  example,  using  the  sample  Cartesian  position 
component  uncertainties  listed  above,  one  can  state:  There  is  a  90%  probability  that  these  components  have 
precisions  belter  than  12, 1.1,  and  0.7  metms  in  X,Y  and  Z  respectively.'  Note  that  this  is  linear  error  for 
each  component  and  not  circular  or  spherical  error. 

Further  details  regarding  siadsiical  statements  and  probabilities  can  be  found  in  AQC  Reference 
PubUcation  No.  28  [1971]. 


9.  Summary 


The  guidelines  and  procedures  presented  here  provide  a  detailed  overview  of  the  STARPREP/GASP 
GPS  geodetic  point  positioning  process.  The  specific  versions  of  STARPREP  and  GASP  which  this 
document  addresses  have  each  been  designated  version  1.0.  These  initial  ’production'  versions  are  based  on 
well  documented,  thoroughly  tested,  easily-mmtained,  portable  FORTRAN  77  source  code. 

As  of  October,  1989,  over  200  STARPREPAjASP  point  positions  have  been  estimated  using 
version  1.0.  If  the  procedures  presented  here  are  followed,  very  few  operational  problems  are  expected. 
Users  of  the  STARPREP/GASP  software  packages  are  encouraged  to  communicate  with  the  authors  of  this 
document  if  any  discrepancies  or  anomalies  are  encountered  when  using  the  software  with  these  guidelines. 

Future  upgrades  of  the  STARPREP/GASP  algorithms  are  anticipated.  These  upgrades  would 
involve  refinements  such  as  the  inclusion  of  a  stochastic  zenith  tropospheric  delay  parameter  in  GASP  and 
the  ability  to  process  new  FICA  configuration  Blocks  from  BEPP/CORE  data.  As  significant  changes, 
enhancements  or  modifications  are  made  to  the  algorithms  or  software,  this  set  of  guidelines  will  be  updated 
as  required. 


Acknowledgments.  The  authors  wish  to  express  gratitude  to  their  past  and  present  colleagues  at 
The  Defense  Mapping  Agency  who  contributed  to  the  development  of  the  procedures,  algorithms  and 
software  describ^  in  this  report.  Most  notably,  these  include  John  Bangert,  Gail  Cherochak,  Barbara 
DeNoyer,  Carol  Finn,  Marilynn  Ison,  Seth  Israel,  Brett  Merritt,  Theodore  Meyer,  Richard  Moore,  Frank 
Mueller,  Maria  Ortiz,  Robert  Pereira,  James  Slater,  Denise  Stutzman,  Brian  Tallman,  George  Tennis  and 
Cesar  Yambot.  Moreover,  the  cooperation  of  coUcagucs  at  the  Naval  Surface  Warfare  Center  is  sincerely 
appreciated. 


30 


UST  OF  REFERENCES  !  BIBUOGRAPHY 

AOC  Reference  Publication  NO.  28,  User's  Guide  to  Understanding  Chart  and  Geodetic  Accuracies,  US  Air 
Force  Aeronautical  Chart  and  Information  Center,  St  Louis,  Mo.,  Sq>teinber  1971. 

Ayrandjian,  A.,  User's  Guide  to  PR7N,  Defense  Mapping  Agency  Systems  Center,  1989. 

ARL,  UT  (Applied  Research  Laboratories,  University  of  Texas  at  Austin),  FIC  File  Document.  Preliminary 
Version,  August  21,  1987. 

Darnell,  A.R.,  User's  Guide  to  GESAR,  Version  1.6,  Naval  Surface  Warfare  Center,  10  November  1987. 

Darnell,  AR.,  Interpolation  /  Propagation  of  Satellite  Clock  Truth  Files,  Naval  Surface  Warfare  Center, 
Internal  Report,  November  1. 1988. 

DMA  Technical  Report,  Department  of  Defense  World  Geodetic  System  1984  -  Its  Definition,  and 
Relationships  with  Local  Geodetic  Systems,  DMATR  83S0.2,  September  30, 1987. 

Evans,  A.G.,  'Comparison  of  GPS  Pseudorange  and  Biased  Doppler  Range  Measurements  to  Demonstrate 
Signal  Multipath  Effects',  Proc.  of  the  Fourth  Intemtl  Geodetic  Symposium  on  Satellite 
Positioning,  Austin,  Tx,  April  28  •  May  2,  1986. 

Fell,  P.J.,  Geodetic  Positioning  Using  a  Global  Positioning  System  of  Satellites,  The  Ohio  State 
University,  Department  of  Geodetic  Science  and  Surveying,  Reprat  Number  299,  Columbus, 
Ohio,  June,  1980. 

Gouldman,  M.W.,  Hermann.  BR..  and  Weedtm,  DI..,  Evaluation  of  GPS  Production  Ephemeris  and  Clock 
Quality.  Proc.  Fifth  International  Symposium  on  Satellite  Positioning.  Las  Cruces,  NM.  March 
1989. 

Hermann,  B.R.(a).  Absolute  Positioning  Results  From  The  Global  Positioning  System,  NSWC  Internal 
Report,  January  1988. 

Hermann.  B.R.(b),  Allan  Variance  Computation  From  OMNIS  Satellite  Clock  Solutions,  NSWC  Internal 
Report.  October  24, 1988. 

Ison,  M.  Preprocessing  for  GPS  Tracking  Data,  DMAHTC  Internal  Memorandum,  1988. 

Malys,  S.<a).  Dispersion  and  Correlation  Among  Transformation  Parameters  Relating  Two  Satellite 
Reference  Frames,  The  Ohio  State  University,  Department  of  Geodetic  Science  and  Surveying, 
Report  Number  392,  Columbus,  Ohio,  July,  1988. 

Malys,  S.(b),  'Similarity  Transformation  Between  NAVSAT  and  GPS  Reference  Frames,  AGU  Chqxnan 
Conference  on  GPS  Measurements  for  Geodynamics’,  R  Lauderdale,  Florida,  September  19-22, 

1988. 

Malys.  S..  and  Ortiz,  MJ.,  Geodetic  Absolute  Positioning  With  Differenced  GPS  Carrier  Beat  Phase  Data' 
P.oc.  Fifth  International  Symposium  on  Satellite  Positioning.  Las  Cruces,  NM,  March  1989. 

Malys,  S.  and  Jensen,  P.A..  Geodetic  Point  Positioning  With  GPS  Carrier  Beat  Phase  Data  From  the 
CASA  UNO  Experiment,  Geophysical  Research  Letters,  CASA  UNO  Issue,  in  print,  1990. 

Rapp,  R.H.,  Geometric  Geodesy.  Part  I,  Dept  of  Geodetic  Science,  The  Ohio  State  University,  Columbus, 
March  1984. 

Swift,  E.,  NSWC's  GPS  Orbit  /  Clock  Determination  System,  Proceedings  of  the  First  International 
Symposium  on  Precise  Positioning  with  the  Global  Positioning  System,  Rockville,  MD,  1985. 

Swift,  E.(a),  Mathematical  Description  of  the  GPS  Multi-Satellite  Filter  /  Smoother,  Naval  Surface 
Warfare  Center,  TR  87-187,  October,  1987. 

Swift,  E.,  'Determination  of  GPS-Realized  WGS  84  Station  Heights  for  the  Combined  Air  FOrce  and  DMA 
Tracking  Networks',  Proc.  Fifth  Intemt'l  Symp.  on  Satellite  Positioning,  Las  Cruces,  NM,  March 

1989. 

Yambot,  C.,  Standard  Operating  Procedures  for  preprocessing  GPS  tracking  Data,  Internal  Memorandum, 
DMAHTC,  1988. 


31 


Appendix  A.  Field  Logs  and  File  Management  Sheet 
GPS  HELD  LOG  INSTRUCTIONS  (OCTOBER  1988) 

This  GPS  Field  Log  will  be  used  with  all  GPS  receivers,  and  serves  as  the  official  record  of 
observations.  The  form  has  been  designed  to  provide  complete  documentation  of  each  tracking 
session,  allowing  the  analyst  to  identify  the  data  collected,  the  operational  scenario  and  any 
notable  events  that  occurred  during  tracking. 

The  form  is  organized  as  follows: 

Tracking  Session  Identification  (top  of  each  page) 

Weather  Data 

Station  Identification  Real-time  Receiver  Position  &  Clock  Solutions 

Operator  Identification  Chronology  of  Events  Log 

Receiver  Information  GPS  Antenna  Height  Worksheet 

Tracking  Scenarios 

Complete  this  form  for  each  tracking  session.  A  "tracking  session*  is  defined  as  a  continuous 
collection  of  data  during  which  there  are  no  intentional,  deliberate  or  planned  breaks  in  the 
tracking.  It  begins  at  a  specific  time  when  recording  begins  and  it  ends  at  a  specific  time  when 
recording  is  terminated.  Any  number  of  changes  in  scenarios  may  occur  during  one  tracking 
session.  Tracking  sessions  differentiate  data  sets  taken  on  the  same  day. 

For  items  that  are  Not  Applicable,  enter  *N/A*. 

For  items  that  are  Unknown,  enter  "Unknown*. 

Enter  page  numbers  at  the  bottom  of  each  sheet 

TRACKING  SESSION  IDENTIFICATION  (ALL  PAGES) 

PROJECT  ID;  Name  of  the  project  or  campaign.  Example:  JPL  Casa  Uno  Campaign  Jan  88. 

STATION  NO.  (DMA):  Station  number  assigned  by  DMA.  If  unassigned,  leave  blank  for  later 
assignment  during  processing.  Example:  85269. 

DATE  (UT):  Day-Month-Year  for  the  Universal  (UT)  GMT  date  when  the  tracking  session 

began. 

Example:  23  Jul  88. 

DAY  of  YEAR:  Day  number  for  the  Universal  (UT)  GMT  date  when  the  tracking  session 
began. 

Examalfl:  209. 

SESSION  NO.:  Letter  code  to  distinguish  one  tracking  session  from  another  on  the  same  date 


1  st  tracking  session .  A 

2nd  tracking  session .  B 

3rd  tracking  session .  C ,  etc. 


NO.  OF  CASSETTES:  Total  number  of  cassettes  used  during  this  tracking  session.  (Cassettes 
should  be  numbered  sequentially  for  a  given  tracking  session,  station  and 
date.) 
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DATA  FILE  NAME; 


File  naming  convention  for  data  stored  in  computer  or  on  diskettes.  Use 
an  1 1  -character  filename  in  the  form  TXXXXYSF.DOO*  where 


T  ■  file  type :  ASCII .  A 

Binary . .  B 


XXXX  a  rightmost  4  characters  of  DMA  station  no. 

Y  >  last  digit  of  the  year 
S  s  tracking  session  no. 

F  <■  disk  file  no.  :  1st  file  this  tracking  session _ 1 

2nd  file  this  tracking  session _ 2 ,  etc. 

DDD  -  Day  no.  (UT)  for  this  session 

Example :  B52698A1.209 


STATION  IDENTIFICATION 

ALTERNATE  STATION  NO.:  Note  other  station  numbers  assigned  to  this  station  (benchmark,  antenna 

site)  either  before  or  during  this  project  Example:  1 1 761 . 

STATION  MARK  AS  STAMPED;Record  the  stamped  information  (not  original  castings)  and  type  of  mark 

(brass,  iron  rod,  cross  chiseled  in  stone,  unmarked,  etc.). 

Example:  Herndon  Optrack  Army  Map  Service  1 965. 

MONUMENT  ESTABLISHED  BY:Organization  responsible  for  establishing  the  monument;  use  the  station 

description  card  if  available;  otherwise  use  the  casted  identification  or 
agency  that  set  the  mark. 

Examrie:  U.S.  Corps  of  Engineers. 

LOCAL  DATUM:  Local  datum  associated  with  tfie  station  mark  and  description  card. 

Example:  NAD  27. 

LOCATION:  Address,  site,  town,  province,  state,  country,  as  appropriate. 

Example:  DMA  Herndon  Facility,  925  Springvale  Rd.,  Great  Falls, 

Fairfax  County,  Virginia,  U.S. A. 

LATfrUDE:  Station  latitude  in  degrees-minutes-seconds;  Indicate  North  (N)  or  South  (S). 

Use  any  available  source  -  map,  receiver  solution,  geodetic  control  sheet,  etc. 

Example:  N  31  deg  25  min  13  sec. 

LONGITUDE:  Station  longitude  in  degrees-minutes-seconds;  Indicate  East  (E)  or  West  (W). 

Use  any  available  source  -  map,  receiver  solution,  geodetic  control  sheet,  etc. 

Example:  E  243  deg  10  min  44  sec. 

HEIGHT ;  Station  height  in  meters.  Enter  either  ellipsoid  height  or  mean  sea  level  height 

as  used  for  receiver  input 

MSL  «  estimated  mean  sea  level  height  (can  use  map) 

ELLIPSOID  -  height  above  the  reference  ellipsoid,  if  known. 

Example:  56.2  m. 
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(XXDRDNATE  SOURCE:  Source  of  latitude,  longitude  and  height  coordinates.  If  map  sheet,  record  series, 
date,  agency,  scale;  If  receiver,  give  manufacturer  and  model. 

Examples:  ONC  E-1 9,  S/76,  DMAAC,  1:1,000,000  or  mission  planning  package. 

REFERENCE  DATUM:  Datum  to  which  the  above  coordinates  are  referred.  Example:  WGS  84. 

IS  THIS  A  BASE  (REFERENCE)  STATION?:  Enter  yes,  no,  or  N/A  as  follows  - 

Yes....  This  station  is  a  primary  survey  point  or  control  point  and  will  be  used  as  one 
of  the  primary  geodetic  positions  for  processing  the  data. 

No .  This  station  is  a  secondary  survey  point  not  used  as  a  primary  reference  station 

for  geodetic  coordinates. 

N/A...  Not  Applicable.  This  station  is  an  independent  point  being  computed  as  an 
absolute  position. 

OTHER  STATIONS  TRACKING  SIMULTANEOUSLY: 

Station  numbers  of  other  stations  simultaneously  tracking  during  this  session. 

ANTENNA  LOCATION:  Antenna  height  in  meters,  as  computed  on  the  GPS  Antenna  Height  Worksheet. 

This  value  is  entered  into  the  receiver.  Antenna  height  is  measured  from  the  top  of  the  station  mark  to 

the  base  of  the  GPS  antenna.  This  does  not  include  the  distance  between  the  electrical  canter  of  the 

antenna  and  the  base.  Example:  1.68  m. 


OPERATOR  IDENTIFJCATION 

OPERATOR'S  NAME:  Name  of  operator  of  receiver.  Example:  I.  M.  Smart. 

OPERATOR'S  ORGANIZATION:  Organization  with  which  the  operator  is  affiliated. 

Example:  DMA. 

RECEf'/ER  INFORMATION 

RECEIVER  MANUFACTURER,  MODEL,  SERIAL  NO.:  Identify  completely  the  receiver  unit  used. 

Example:  Texas  Instruments,  Tl  41 00,  SN  850005551 83. 

ANTENNA  MANUFACTURER,  MODEL,  SERIAL  NO.:  Identify  completely  the  antenna  used. 

Example:  Texas  Instruments,  TI4100,  SN4141. 

CLOCK  MANUFACTURER.  MODEL,  TYPE,  SERIAL  NO.:  Identify  completely  the  external  oscillator 

used,  if  applicable. 

Example:  Efratom,  FRK-30,  Rubidium,  SN  8585. 

_ ^NO  E)CrERNAL  OSCILLATOR  USED: 

If  no  external  oscillator  was  used,  check  here  on  the  form. 

NAVIGATION  PROCESSOR  SOFTWARE  NAME  AND  VERSION: 

Name  and  version  of  navigation  processor  software  used 
inside  the  receiver  to  track  satellites  and  make 
measurements. 

Example:  GESAR  Version  1 .9. 
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MEASUREMENT  INTERVAL 


MEASUREMENT  RECORDING  INTERVAL 

Rate  (interval)  at  which  measurements  are  recorded  or  written 

to  tape  or  disk.  Example:  30  sec. 

REAL-TIME  SOLUTION  INTERVAL 

Rate  (interval)  at  which  the  real-time  navigation  solution 

(position)  is  computed  or  updated  by  the  receiver. 
Example:  30  sec. 


TRACKING  SCENARIOS 

For  each  scenario  provided  by  the  lead  organization,  enter  the  following  information: 

START  TIME  (UT)  :  Enter  the  SCHEDULED  start  time  and  the  ACTUAL  start  time. 
Example:  1020  1030. 

SATELLITES  TRACKED:  Enter  the  PRN  numbers  of  the  satellites  tracked. 

Example:  6  9  11  13. 


WEATHER  DATA 

At  planned,  regular  intervals,  record  the  following  weather  data  and  key  into  the  receiver  if 
required: 

TIME  (UT) :  UT  time  of  the  weather  station  readings.  Example:  1430. 

TEMPERATURE  (°C) :  Wet  bulb  and  dry  bulb  temperatures  in  degrees  Celsius. 

Example:  24.6/36.8. 

PRESSURE  (Mb) :  Barometric  pressure  in  millibars.  Example:  984. 

RELATIVE  HUMIDITY  (%)  :Relative  humidity  in  percent.  Example:  59. 

REAL-TIME  RECEIVER  POSITION  &  CLOCK  SOLUTIONS 

During  the  tracking  session,  it  is  useful  to  record  the  position  and  clock  solutions  in  order  to 
monitor  changes  in  the  receiver's  position  and  dock  estimates  for  quality  control  purposes. 

This  should  be  done  at  planned,  regular  intervals. 

TIME  (UT)  :UT  time  at  which  the  real-time  receiver  position  and  clock  solutions  are  entered  on 
the  form.  Example:  2330. 

LATITUDE  (deg,  min,  sec)  ;Latitude  displayed  by  the  receiver  in  degrees,  minutes  and  seconds. 
If  the  receiver  only  displays  positions  in  decimal  degrees,  use  decimal  degrees. 

Example:  62  -  25  -  40  or  62.4278. 

LONGITUDE  (deg,  min,  sec):Longitude  displayed  by  the  receiver  in  degrees,  minutes  and 
seconds.  If  the  receiver  only  displays  positions  in  decimal  degrees,  use  decimal  degrees. 

Ex^ple:  288-30-18  or  288.5050. 
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HEIGHT  (m)  :Ellipsoid  height  displayed  by  the  receiver  in  meters.  Example:  1 41 .3 


CLOCK  BIAS  (time  units) :  Receiver  clock  offset  from  GPS  time  displayed  by  the  receiver  in 
seconds.  Example:  -0.1407  microsec. 

CLOCK  DRIFT  (time/time  units)  :Receiver  clock  drift  (frequency  offset)  displayed  by  the 
receiver  in  sec/sec,  microsec/hr,  etc. 

Example:  0.7671  microsec/hour. 


CHRONOLOGY  OF  EVENTS  LOG 

TIME  (UT)  :UT  time  at  which  an  event  occurred.  Example:  1 31 7. 

COMMENTS:  Record  error  messages.  Describe  special  actions  taken,  note  if  and  when  the 
receiver  was  left  unattended,  document  interruptions  to  tracking,  etc. 


GPS  ANTENNA  HEIGHT  WORKSHEET 

Compute  the  height  of  the  antenna  base  above  the  mark  as  indicated  on  the  worksheet  Measure 
heights  in  both  meters  and  inches.  Use  the  antenna  constants  provided  on  the  worksheet  or 
supply  new  ones  as  appropriate  for  the  antenna  being  used.  Record  any  new  constants  on  this 
worksheet.  Transfer  the  ‘Computed  Antenna  Height*  to  page  1  of  this  log  form  under  the  item 
labeled  ‘Antenna  Location*. 
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GPS  FIELD  LOG 


Station  No.  (DMA) 

Date{UT) _ 

Day  Number _ 

Session  No. _ 

No.  of  Cassettes_ 

Project  ID _  Data  File  Name _ 

STATION  IDENTIHCATION 


Alternate  Station  Number _  Latitude _ deg _ ^min _ ^sec 

Longitude _ deg _ min _ ^sec 

Station  Mark  as  Stamped _  Height  (meters) 

_  MSL _ 

Monument  Established  by _  Ellipsoid  _ _ 

Local  Datum _  Coordinate  Source _ 

Location _  Reference  Datum _ 


Is  this  a  base  (reference)  station?  _ Yes  _ No  _ N/A 

Other  stations  tracking  simultaneousi}': 

Station  No. _  Station  No. _  Station  No.. 

Station  No. _  Station  No. _  Station  No.. 

Antenna  Location 

Antenna  Height  (meters) _ 

(see  attached  worksheet) 

OPERATOR  IDENTIFICATION 

Operator's  Name _  Operator's  Organization _ 


RECEIVER  INFORMATION 

Manufacturer  Model  Type  Serial  No. 

Receiver  _  _  XXXJCXXXXX  _ 

Antenna  _  _  XXXXXXXXX  _ 

Clock  _  _  _  _ 

_ No  external  oscillator  used 
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10/88 

GPS  FIELD  LOG 


Project  ID _ 

Navigation  Processor  Software 

Name _ 

Measurement  Interval 
Measurement  Recording  Interval 

TRACKING  SCENARIOS 


Page _ ^of 

Station  No.  (DMA) _ 

Date  (UT) _ 

Day  Number _ 

Session  No. _ 

No.  of  Cassettes _ 

Data  RIe  Name _ 


Version 


Real-time  Solution  Interval 


Start  Time  (UT) 

Scenario  Scheduled  Actual 

1  _  _ 

2  _  _ 

3  _  _ 

4  _  _ 

5  _  _ 

6  _  _ 

7  _  _ 

8  _  _ 

9  _  _ 

10  _  _ 


Satellites  Tracked 
PRN  Nos. 


WEATHER  DATA 


Time  (UT) 


Temperature  (®C) 

Wet  Bulb  Dry  Bulb 


Relative 

Pressure  (Mb)  Humidity  (%) 


Page _ of 


38 


GPS  FIELD  LOG 


Project  ID_ 


Station  No.  (DMA) 

Date(UT) _ 

Day  Number _ 

Session  No. _ 

No.  of  Cassettes_ 
Data  Rie  Name 


REAL-TIME  RECEIVER  POSITION  &  CLOCK  SOLUTIONS 

Time  Latitude  Lonptude  Height  Clk  Bias  Clk  Drift 

(UT)  (deg,inin^ec)  (deg,min^)  (m)  (sec)  (sec/sec) 


CHRONOLOGY  OF  EVENTS  LOG 


Time  (UT) 


Comments 


Station  No.  (DMA). 

Date  (UT) _ 

Day  Nuinber _ 

Session  No. _ 

No.  of  Cassettes 
Data  File  Name _ 


Orientation  to  North 


BASE  OF  PREAMP 


ComputedAntenna 

Height— 

A  = _ meters 

_ inches 


A  =  (C 


Ht.  of  Mark 
aboveGround 

D  = _ meters 

_ inches 


Page _ of 
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Appendix  B.  Sample  STARPREP  Output 


The  selected  pages  of  STARPREP  output  reproduced  here  represent  the  typical 
output  from  a  normal  program  execution.  The  entire  output  has  not  been  included  since 
much  of  it  is  simply  repeated  for  each  satellite  which  was  processed.  The  complete  output 
from  which  these  samples  were  extracted  was  97  pages. 
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o  o  0990  o 

4  •  «  4  •  • 

Ml  Ok  Ok  Ok  Ok  »  > 

IM  ^099^0 

C  va  ookOkOkO^ 

2  E  cS&fe&ftte 


imo  ovjQOao 

Ml  o  o  o  o  o  o 


«•  A«  n  O  kA  4 


fk«  <4  -4  flO  kA 
Ok  4»  lA  M*  ^  O 
^4  rj  -'aiA  -4 


^1  4  Ok  ••  A|  *A 


O  O  k^  kA 


a. 


««UltCL»SSJf  1  5-**  ilAkJ'LF  rL»nVli.Ll. 


vs  Tint 


1 


SATULllE  J 

.U270-0U2 

•fcoct-oa; 

.7e5e-oo2 

.76S2-007 

.7<«S-OU2 

.7240-C02 

.7CJS-OJ2 

• 68fc8*0U2 

•  (i62Z-OU2 
.t)416-C02 

•  6210-0U2 
.faOU4-0(<2 
.5798-002 
.5592-002 
.5386-002 
.5181-002 
.4975-002 
.<»7to9-ou2 


NE1.A1  IV11T 

I - , - - - 


(IXXX 


XXX 


XX 


XXX 


XX 


I 

I 

1  XX 

-  XX 


XX 

X 


I  XX 

1  X 

-  » 


I  X 

I  X 

I  I 

1  XX 
-  X 


xxxxxxxxxxxxx 

xxxx 

XXX 


.SU4S«30t  .5C67«U)6  .S089AU06  .51  11 4C0b  .51334006 


**UHCLASSXri£0**  STARfREP  PLATTVlLLC  DAT  223  1929 


SATELLITE  3 

LARTH  KOlAtE 

vs  TIHL 

-.8979-002 

-.1005-001 

XX 

-.1111-001 

X 

-.1218-001 

XX 

-.1325-001 

- 

XX 

-.1< 32-001 

XX 

-.1538-001 

XX 

-. 1643-001 

XXX 

-.1752-001 

XX 

-. 1859-001 

• 

XXX 

-.1966-001 

X 

-.2072-001 

XXX 

-.2179-001 

• 

xxx 

-.2286-001 

XXX 

-.2393-001 

xxxx 

-.2499-001 

xxxx 

-.2606-001 

-.2713-001 

]  Axxxnx 

•  XXXIXIXX 

.50454006  .50674006 

.S0894W6  .51  114006  .51334006 

••UNClASSirX  E0«*  S1ARFRLP  PLATTVILLE  OAT  223  1929 


satellite  3 

SAI. 

M  1 ENNA 

vs  UBE 

.8759-003 

XX 

.6751-003 

1 

xxx 

.8742-003 

1 

xxx 

.8734-003 

1 

xxx 

.6726-003 

• 

xxx 

.6716-003 

1 

xxx 

.8709-003 

) 

xxx 

.6701-003 

1 

A  XX 

.6693-003 

1 

xxx 

.8664-003 

- 

xxx 

.8676-003 

I 

xxx 

.8668-003 

I 

xxxx 

.8660-003 

- 

xxx 

.8651-003 

1 

xxx 

.8643-003 

1 

XAX  A 

.8635-003 

1 

xkxx 

.8627-003 

1 

VX  XX 

.8618-003 

-  %%% 

.50454006 

.S0674U06 

.50894006  .51  11  4006  .St33*(i06 

5  5 


SATeiLZTC  i  NOPFLt  TROfO  VS  TiHC 

.9577-002  -  t 
.9167-002  1  X 
.6760-002  1  X 
.8354-002  1  » 

.7948-002  -  X 

.7541-002  1  X 

.7135-002  1  X 

.6728-002  I  XX 

.63  2  2  -002  I  XX 

.5916-002  -  XX 

.5509-002  1  a 

.5103-002  I  XX 

.4697-302  -  XXXX 

.4290-002  I  XXXX 

.3884-002  1  XXXXX 

.34  77-002  1  XXXXXXX 

.3071-002  I  XXXXXXXXXX 

.2665-002  -  X 


.4657-002  -  XXXX 

.4290-002  I  XXXX 

.3884-002  1  XXXXX 

.34  77-002  1  XXXXXXX 

.3071-002  I  XXXXXXXXXX 

.2665-002  -  XXXXXXX 

.5045*006  .5067*006  .5U89*U06  .5111*006  .5133*006 

SATCLL11C  6  HAS  t.C  PIT  TlNC  WITHIN  THE  TOLERANCE  TOR  100  MLASUREHENTS 

bETHECt*  WEEK  SCC  TlHl  49 79 09.013u66l  26  An6  WEEK  500  TINE  5  0  0969.010  3324  42 

SATELLITE  6  MAS  kC  PET  TINE  WITHIN  THE  TOLERANCE  TOR  100  NLASURCNEUTS 

BETWCEN  WEEK  5CC  TIHC  SCQ999.0102860U9  AN6  WEEK  SOU  TINE  50 3969 . o04 1 8042 1 

satellite  6  HAS  KC  PET  TINE  WITHIN  THE  TOLLRANC  E  FOR  IS  HEASURCtlCNTS 

UCTtfECN  WELK  SLC  TINE  503999. 004106407  AN6  WEEK  SOU  TINE  SO 4359. U03 204602 


•*UNCLASS1T1CD**  STARTRLP  PLATT91LLE  OAT  223  1989 


SATELLITE  6  RANGE  lOKO  V$  TINE  _ _ ^ 

-.2295-002  -  X  X 

-.2766-002  I  *  X 

-.3237-002  I  * 

-.3708-002  I  X 

-.4179-002  -  X  XXX 

-.4650-002  I  I  XX  XX  X  X 

-.5121-002  IX  XXX  XXX 

-.5592-002  IXX  XXXXX  XXXX 

-.6063-002  1  XX  X  K  XX 

-.6534-002  -  *  X  X  XX  X  X 

-.7004-002  I  XXX  X 

-.7475-002  1  X  XXX 

-.7946-002  -  XXX 

-.4417-002  I  X  X  XXX 

-.8888-002  1  X 

-.9359-002  I  XXX 

-.9830-002  I 

-.1030-001  -  _  _  *  j 

.4980*006  .4996*006  .5012*006  .5028*006  .5044*006 


)£LL11E  6  UA1.  (LOCK  VS  Tint 

.10184005  -  KX  *  * 

.101840113  1  XXI 

>  10184003  1  XXX 

.10184003  I  XXX 

.10184003  •  xxx 

.10184003  I  xxx 

.10184003  I  XXX 

.10184003  1  XXX 

.10184003  I  XXX 

.1018  4003  -  xxxx  I 

.10184003  I  AXXX 

.10584003  1  XXX 

.10184003  -  'XXIX 

.10184003  I  xxx 

.10584003  t  XXX 

.1058*003  1  xxxx 

.10184003  I  xxxx 

.10184003  -  I 

.X98040pe  .«99(i4X6  .30124006  .5u2s4;;C6  .S0X44u0b 


•*UHCL*SSiri  E0*4  ITARfXLP  PLXTTV1LL8  DAY  223  19(9 


SAT  LLlITE 

13A6-C 

-.1A93-C 

-.16A0-C 

1787-t 

-.1933-1 

-.20b0-C 

-.i227-C 

-.2374-C 

-.2S21-C 

-.2668-( 

-.2815-t 

-.i962-C 

-.31C8-C 

-.3255-( 

-.34o2-t 

-.3S49-{ 

-.1696-C 

-.38a3-( 


VS  Tint 


6  RCLATlVIir  VS 

1 - 1 - 

2  -  XX 
2  1  XXI 
2  1  XX 

2  1  xxx 

2  -  XX 

2  I  xxx 

2  I  XXXX 

2  1  xxx 


02  1  xxx 

02  1  xxx 

02  '  xxx 

02  1  X  XX  X 

02  1  xxxx 

02  1  xxxx 

02  1  xxxxx 

(l2  *  xxxx 

.49804008  .49964006  .50124006  .10284006  .50444006 


••UhCLA^Si  F  1  ta  ••  :.UI>FN|.P 


PLArrviLH.  o*»  Hi  I9^v 


SATtLLlTE  6  SAT.  ANTLHItA  VS  Tint 

.eSS7-003  -  SXAU 

.b84S-0b3  I  XMXXA 

.8S13''003  1  XXXX 

.8811-003  1  XXXX 

.88C9-OU3  -  XXX 

.8797-003  I  XX 

.8783-003  I  XXX 

.8773-003  I  XX 

.8761-003  1  IXX 

.8750-003  -  XX 

.8738-003  I  XX 

.87.6-003  1  XXX 

.87lt-0U3  -  XXXX 

.8707-003  I  XXX 

.8690-003  I  XXXX 

.867t-Cu3  I  XXX 

.8666-003  1  XXX 

.«65*-003  -  XX 

.498O*00l  .499o*C06  .5012«006  .50784006  .5044*006 


L 


••UNCLASSiri  10**  SlAXXRLP  PLATTVltLE  OAT  773  1979 


c 


c 

c 

c 


SATLttlTE  6  HOPFLC  TAOPO  VS  Tint 


.4637-007 
.4483-007 
.4334-007 
.4185-007 
.40  3  6  -00  7 
.3887-007 
.3738-007 
.3589-007 
.3439-007 
.3290-007 
.3141-002 
.7992-007 
.7843-002 
.7634-002 
.2545-002 
.2396-007 
.2247-002 
.2098-002 


1 


X 

X 

XX 

X 

X 

XX 

X 

XX 

X 

XX 

XXX 
XXXX 
XX  XX 
XXX 

xxxxx 

xxxxxx 

xxxxxxxxx 

XXXM  xxxxx 


.4980*006  .4996*00  6  .5017*^6  .50  78*006  .5044*006 


SATELLITE  9  HAS  AC 
eETwEEN  WEEK  5CC 


PIT  Tint  WITHIN  THE  TOIEXANCE  FOX  100  HEASUXEnENTS 

TJHC  49  79  09.008692891  ANX  LEEK  500  TINE  5U6309. 008059484 


SATELLITE  9  HAS  AC 
HETHEEN  WEEK  SCC 


PIT  Tint  WITHIN  THE  TOLEXANCE  EOX  iQO  HtASUREnEHTS 

Tint  506339.008001427  ANX  LEEK  500  TIHe  5o9339. 001 020872 


satellite  9  HAS  AC 
OETWEEN  week  5CC 


PET  T1F«E  WITHIN  THE  T0LE*ANC  E  FOX  54  HEASUXEnENTS 

TINE  5  0  93  69.00094  3164  AAX  WEEK  500  TINE  510958.996  764535 


O 


••UACLASSIEIEO**  S1AXFXEP  PLATTVILEE  XAT  723  19E9 


L. 


SATELLITE  9  RANGE  ICAO  VS 

-.2738-007  -  X 
-.34  0  7-002  1  X 

-.4076-002  1 
-.4746-007  I  X  X 
-.5415-002  -  I 
-.6084-002  I  X 
-.6753-002  I  X 

-.7477-002  I  X 

-.3092-002  1  X 
-.8761-002  - 
-.9430-002  I 
-.1010-001  I 
-.1027-061  - 
-.1144-001  I 
-•l2l1-00l  I 
-.1277-001  I 
-.1343-Oul  I 
-.1411-001  - 

1 - j - 

.49Eu*U06  .5012*606 


INE 

X 

I  X 

8  A 

X  X  XX  XX  X  X 
XXX  X  X 
XXX 

.  X 

XX  X  XX  X 
AX  X 

X  XX  XX  X 
X  XXX  X 

X 

XX 

X 

.5044*006  .5077*006  .5109*006 


X 

XX 

X 


# 


5  8 


**UUCLASS1  FI  EU**  STARPRLP 

PLATTVILLE 

DAY  223  1989 

DATA  CORRECTIONS  LtMMARY 

PRN  NO. 

NO.  OES.  INPUT 

PILE  NO. 

INPUT  FILE  NAME 

3 

247 

11 

Gt0STAP*05  35FLA9  223T. 

0 

163 

12 

GE0STAR*06S5FLA9  223T. 

9 

EGA 

13 

GEOSTAR*Gy  25PLA9  223T. 

1  1 

•<.62 

14 

G£OSTAfi*n  85FlA9223T. 

12 

AlS 

15 

GEOST  aR*12  E5FLA9  223T. 

13 

243 

16 

6£0STAR*13  25  FuA9223T. 

EPHtHLRIS  SOURCL:  RLFERLNCE 


PRN  NO. 

tPHEF.  FILL  NO. 

EPHEK.  FILE  NAME 

3 

21 

GPS*EF03ii921b. 

6 

2  2 

&PS*EF06ij921  &. 

0 

-  •» 

GPS*EFU9 *921£  . 

11 

24 

GPS*EF11S9216. 

12 

^  5 

GPS*EF12i921  a. 

13 

26 

GPS*£F13a92l6. 

CORRECTIONS  : 

CCKPLTED 

APPLIED 

RANGE  lONO 

range  iono 

DCPPLER  lONO 

DOPPLER  ICNO 

SAT  .  CLOCK 

SAT.  CLOCK 

R  EL  ATIVITY 

RELATIVITY 

earth  rotate 

EARTH  ROTATE 

S  AT  .  antenna 

SAT.  ANTENNA 

CHAC  TROPO 
hCPFLO  TROPO 

HOPFL&  TROPO 

OUTPUT  FI4.E  NO.;  1C 

OUTPUT  file  NAME:  GL  C 'T  AR*P  TP  LA9223P  . 


♦•UJ.CEAESi  FI  Ei**  ETAFFRlP  PLATTVIELw  T  AY  223  195-/ 


M  M  r*  K> 
rJ  I'j  rj  rt 

IM  IM  (\J  l>J 


••u»,CL*4Sif  1  tj** 


rLATtVlLH  DAY  cii  JYff 


rili  SUHKARV 


thc  roiLOuluG  FiLcs  kim  Acasscp  pu»iuc  ihu  iiu 

NO.  IlHL 


E  NO. 

FILE  NAPl 

A 

UEOSIAR  •SIFTlOiauC- . 

7 

GC0$1AR««.A5PLAA9223. 

10 

GE0STA6  <PTFLA9l2JP. 

11 

6COSTAR*OII3rLA922  3T, 

12 

CEUSTAk«l:EE5PLA9223T. 

13 

GC05TAC'a9ESPLA92?3T. 

U 

bE0STAR«1  1i5PLA9223T, 

15 

bCUST  AS  M  2ESPLA9  223T. 

16 

GEOSIA6«13e5PLA9223T, 

21 

OFS*EFC;592U. 

.2 

GFS*EFLCES218. 

23 

GPS*£rC16f2ie. 

2A 

bPS*EF1189212. 

23 

GFS*EF12tiS:ie. 

26 

bFS*EF1!L921H. 

31 

GEOSTAR  •nTESFLA9223P 

32 

GEOST  as  'U  CESPLA9  22  3P 

33 

GEOSTAR  «i  TEIPLA9223P 

36 

GCUSIAR<S1ESPLA9223P 

60 

GPS*PCei21E. 

41 

GEOST  AR<PTt3223La;T. 

42 

GE0STAR*PTE6;23L0UT. 

53 

GEOSTAR«P119223LOUT. 

56 

GE0STAR«PT11223L0UT. 

55 

C  COSTA  (•FT  122231  OUT. 

'•UNCLASSintU**  SJAAfULP  PLATTVKLS  OAT  223 


•*UNCLASSiriEO*»  SURFREP  PLAITVlLLL  DAT  223 
SA  6(UStA««P113;'23LCUT. 


CPEMCD  PURPOSE 


919 


919 


SAILLLITE 

NO. 

3 

satellite 

NO. 

6 

satellite 

NO. 

9 

satell  he 

NO. 

11 

SATELL  ITE 

NO. 

12 

satellite 

NO. 

13 

riCA-rORHATTLD  INPUT  DATA 
POINT  riUE  (CORRLCTEb  NEASURE"ENIS) 
LOiT/TTAO  NEASUREniNISt  SATELLITE  HO. 

EDlT/TTAb  HEASUREPLNTS,  SATCLLETL  MO. 

LOlT/TTAfi  HEASUPEPENTS.  SATELLITE  NO. 

EDIT/TTAO  HEASURLPENTSt  SATELLITE  HO.  11 

ED1T/T1AO  HEASUREPLNIS,  SATELLITL  NO.  12 

EPET/TTAG  HEASUREHEHTS.  SATELLITE  HO.  13 

REFERENCE  EPHEHERIS,  SA1 
REFEREHCL  EPHEHERlSt  SATi 
REFERCHCE  EPHEHERIS.  SATl 
REFCRENtE  EPHEHERIS.  SAT 
REFEREHCE  EPHEHERIS.  SATl 
REFERCHCE  EPHEHERIS.  SAT; 

TRANSACTION  HEATHER  PAIR 
TRAHSACTIOH  HCEEIVER  DATA 
transaction  satellitc  data 
TRANSALTION  STATION  DATA 


PAI 


DAT 


c 

c 

c 

c 

to 


e  0 


This  PhOGHAH  USLS  IHt  POINT  f  Itt  OT  TMC  FBEPROClSSOR 

PkOCDAfl  STkSPMfP  1C  PtiOOulC  *  pOiMT  PoSlTiOk  fOOP  tllf^liCOCLD 
6PS  CARKU*  bt*T  t  PASk  lAlA  (COLl,Lt1lO  WITH  A  TlAlCvi  htUl\|LR) 


o  •  ac  «  oc 


■■•i 


a  o  o  a  o 

o  o  a  o  «A 

•  *  «  •  • 

O  >4  O 


•  ‘Vr'ir'.V/  •. 


X  ^  N| 
W  ^  «  «K 


X  M  ^ 

I  OB  O  W  <  < 

U  U  3  ««  >4 

M  Ik  9  9 

:  ^  mt  a 

A  Id  V  o  o 


*  *A  jt  » 

do  «  < 

:  £  ^  m  A  ^ 

O  X  z 

I  >1  u<  X  <  < 

I  X  ^  ^ 

I  X  a»  X  (/» k/t 


<  ^  su 
>i  «  3 


^  n  4C  « 

«  o 

M  C  X  & 

3  X  O  d 

O  O  M 

UJ  «k  ««  W 

i«i  X  «  a 

W  X 

a  &  *«  .J 

•«i  ,  ^  ^ 

►•  O  v»  •• 

SI  *•  J  3 

•J  o  ^  w 

IW  S4  <  > 

W»  W9  •«  3 

a  »-  a 

u  X  < 

a  w  u/ 

a  3  X 

d  ^33 


«<  «l  X  >* 

a  J  ►- 


W  «  u*  o 

X  A  ^  X 

M  X  W  X 

X  <  >•  a 

IW  ^  w  X 

kte  »/|  < 


o  i^r» 

MO  O  ^ 

C  •  fw  ^ 

O  O  M  O’! 

♦  rw 

o  o  •  ♦  • 

o  •  •-  •• 


M  X 

M  a 

O  fM  »- 

<  MX 


lA  W^  O  M  1^  O^ 

►>  M  •>  a  3 


I  O  ^  ^  M  »A 


o  o  o 

y»  *•  fk*  N* 

x»*  ^ 

X  -A  K» 

^  -«o  o 

3  a«A  ^ 


Ml  ^  O 

53  -.8:3 

X  !•>  43  C^  O 

>^1  3 »  ^ 

ac  KO 
X  ^  ^  ^ 

O  ^  lA  n  ^ 

.1  • 

S«  «  O  V>  ITk 


->  o 
^  ^  o 

X  »  -5 
«l  v>  w 

Ok  • 

*«  W  iW  (X  c 

X  X  >  4 
O  SO' 

«N  Ml  >>  «• 

Oc  O  9^ 

C.  3  •»  • 

«  UI  <C  C 
3  .J  'f 

O 


THE  «L*N  0-C  EOM  7i7  GASP  OU  SEHV AfitE  G  IS  I  -13.075  CH 


•u»ti  or  DLLl*  PiCLroRiiMbE  -  DELIA  CARHIlK  BtAl  PIIASL  fOR  LACH  SATLLLIIE  IN  IHJS  DAlA  Ml 


£  C  C  r  C  iC  . 


:s  •£££££££: 


:££.C£9;£Z:C£££  C£. 


^  O  ^  J'^O*  J*  A*>kj  *^0 

wn  **  r\4  O  «*.  tnorw  a  *•»»•  «v«  O  O  dt*  O  OC  Opo  •#'«rOfw*»*»  ♦it-mCC 


fSM  *n  «•  nO  W'« 


£  ££££££££] 


^CS£SCKXCCCXCS£C££i:CZSt;CC££ 


eaccoowNaofMco*^OM^*r(v*vNtf>iA  0'000'OAO«o<><<o«*^^o^'0<Njrvjfsio<<#'^ 


Ai  S£'S£££3X£S«SSSXXX££SSXSSSZS£SSZS£SXS£££ 


<  «  <  <  <  <  «  < 


r>*  r^«^4  <0  o»  'O 
-3  4  *1 

o<^  ^ 

•4r  ^  o  aa 

fN»  cri  «u  r' 

^  -N*  f-  3  ••  3 


i^-n-nx^^fn  oaA'Q  oa  ai 


tA«AA  aA  aaA  a  ^  a  a  a  a  a  . 
I  I  I  I  <  f  I  I  I  I  <  I  I  I 


^OOOidOOOOOOtf^OOV 

;<<<<<«<«««««««« 

(JSXxxaatacQsxsxxaKS^ 

loooowoucioeeococ 

^aaaaaaaaaaa^aa: 

ssaaa  aa  asaa  aaaa: 

*it%tAy%i/iv>yttn¥%y%t/tv»v%v*y»\ 

DXx&a»&.xa»;kdAX^a>a.^fi 


,  <  <  <  <  ^  X  <  t 


I  1  I  I  I  I  I  I 
j  ^  jj  ^ 

!2^a2£Z££ 

|X£XX£.«£X 

>oQeouc^w(-* 

ssasaa  aaa 


aoAAAAaAa 


J  o  ' 

«  A  ■ 

^  m^  r" 

»A  <0  »0  ^  Ol 

^  ••iW 

V) 


aoaatija  3-3  a3*“3^'3  a-^  ^33 

c.<oowwttowowoo(aooc9&»e^c)ouoooeoob>wu>oukiuji>.>uu2(a* 

•  «••••••••»••••••••••••****•••••«•***** 

^OAOAON^<K<^OAO>0•C^O>OA9><^9‘OAOA<r^<^(^<^CAC^O‘C^^^C^i^O^^C^^O«^OA<^ 

••••••OOO  ••«‘Jmeo  •^iNi|wfMK«weOOnf^#N*^'r»»*OOOw>K  OO 

a  194><,^>JkOD  ^o 

^9A»VN0A0A0A0A|^UCPOO3OOaUOOOOO0OOUUOOOO^«*^^^*>** 

aa3a3o*»aa**f»R»3^^ao  3*»*»i^^^aa«a  aa«»  3*»3aj**^o 

oocoooooooooooooooooooooooooooocoooccco 

0>0‘^0«0^Qa0a^0»0.^O-0a0a0R9>9>Oa^0aC>0>^9^^^^O'CS9a9‘Oa9-^C‘0‘O‘O'^ 

^  ^ 't  >r  Af'i  A'-f  O  30-403  »  30‘J-‘'J 
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Appendix  D.  Troubleshooting 


The  chart  provided  below  is  intended  to  help  new  STARPREPA3ASP  users  diagnose  and  solve  processing 
problems  which  may  occur  during  familiarization  and  training  phases. 


Symptom 

Possible  Problem 

Solution 

STARPREP  defaults  to 
broadcast  ephemendes 

One  or  more  precise 
ephemeris  files  missing 

Include  all  required  ^hemeris 
files  in  runstream 

Tracking  span  is  to  large 
(requiring  over  100  ephemeris  points) 

Shorten  Tracking  span  or 
eliminate  large  ^ta  gap 

Point  file  contains 
meaningless  interpolated 
ephemerides 

or 

GASP  message:  no  data  in 
point  file 

Data  spans  week  crossover, 
but  standard  (full  week)  ephemerides 
were  assigned  in  the  STARPREP 
nmstream 

Use  half-week' 
ephemeris  /  clock  files 

or 

Edit  FICA  file  to  remove 
week  crossover 

Non  specific  enor  message: 
when  running  STARPREP 
on  UNISYS:  Guard  mode 

Wrong  Precise  clock  file 
assigned 

Assign  correct 
clock  file 

GASP  position  plots  do  not 
converge  or 

High  conelation  coefficients 

High  standard  deviations 

Data  set  is  too  short  (or  too  much  rejected) 

or 

Very  poor  geometry 

Use  different  data  set 
or 

Determine  why  too  much  data 
was  rejected  and  lesolve 

GASP  residual  plots  do  not 
appear  random  at  last  iteration 
or 

Large  percentage  (>30%)  of  data 

rejected 

or 

Most  data  rejected  in  GASP  by 
RMS  screening 

A  priori  station  coordinates 
are  too  gross 
or 

Data  set  contains  an  eiratic  satellite 
(check  delta  pseudorange  •  delta 
carrier  beat  phase  biases,  and  rejected 
data  by  sateUiie) 

Rerun  STARPREP  and  GASP 
with  improved  a  priori 
station  coordinates  (obtained 
from  initial  GASP  result  if 
some  data  was  processed) 
or 

Eliminate  erratic  satellite 

GASP  position  plots  do  not 
converge  or 

High  RMS  of  Residuals 

Missing  or  empty  clock  file 

Assign  correct  clock  file 
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